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A recently introduced method utilizing dimensional continuation is em- 
ployed to compute the energy loss rate for a non-relativistic particle moving 
through a highly ionized plasma. No restriction is made on the charge, mass, 
or speed of this particle. It is, however, assumed that the plasma is not 
strongly coupled in the sense that the dimensionless plasma coupling parame- 
ter g = e^Ko/^T^T is small, where Ko is the Debye wave number of the plasma. 
To leading and next-to- leading order in this coupling, dE/dx is of the generic 
form g"^ ln[C(7^]. The precise numerical coefficient out in front of the loga- 
rithm is well known. We compute the constant C under the logarithm exactly 
for arbitrary particle speeds. Our exact results differ from approximations 
given in the literature. The differences are in the range of 20% for cases rel- 
evant to inertial confinement fusion experiments. The same method is also 
employed to compute the rate of momentum loss for a projectile moving in 
a plasma, and the rate at which two plasmas at different temperatures come 
into thermal equilibrium. Again these calculations are done precisely to the 
order given above. The loss rates of energy and momentum uniquely define 
a Fokker-Planck equation that describes particle motion in the plasma. The 
coefficients determined in this way are thus well-defined, contain no arbitrary 
parameters or cutoffs, and are accurate to the order described. This Fokker- 
Planck equation describes the straggling — the spreading in the longitudinal 
position of a group of particles with a common initial velocity and position — 
and the transverse diffusion of a beam of particles. It should be emphasized 
that our work does not involve a model, but rather it is a precisely defined 
evaluation of the leading terms in a well-defined perturbation theory. 
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I. INTRODUCTION 



The methods of quantum field theory (QFT), originally developed to describe the inter- 
actions of elementary particles, have since been successfully applied in several other fields 
of physics, including plasma theory. Condensed matter physics in particular has been a 
proving ground for the utility of QFT methodology. The Kondo problem, phase transfor- 
mations in magnetic metals, and a host of other problems in condensed matter theory have 
been tackled using QFT methods [1-3]. Classical non-equilibrium reaction-diffusion sys- 
tems have been studied using QFT techniques since the pioneering work of Doi in 1976 [4]. 
More recently the methods of effective field theory were applied to multicomponent, fully 
ionized plasmas by Brown and Yaffe [5]. These methods significantly simplify high-order 
perturbative calculations and clarify the structure of the theory. Here we use the method of 
dimensional continuation, which was originally developed as a means of regularizing diver- 
gent integrals that arise in perturbative calculations in gauge-invariant QFT, to calculate the 
stopping power and temperature equilibration in highly ionized plasmas. The well-known 
Lenard-Balescu kinetic equation [6,7] describes the long-distance, collective excitations of 
the plasma, whereas the Boltzmann equation for pure Coulomb scattering describes the 
short-distance, hard collisions of the plasma particles. A complete description of the plasma 
includes both the long- and short-distance physics encoded in the Lenard-Balescu and Boltz- 
mann equations, but each contribution is divergent if integrated over all of three-dimensional 
space. A finite sum, the large Coulomb logarithm and its coefficient, can be obtained by 
introducing cutoffs, but this approach docs not determine the constants that accompany the 
logarithm. These additional constants arc given by the convergent kinetic equation method 
[8], but spurious higher-order terms arc introduced. In contrast, our method, which is based 
on a rigorous expansion in a small parameter and dimensional continuation of the Lenard- 
Balescu and Boltzmann equations, gives the constants accompanying the Coulomb logarithm 
but no spurious higher-order terms. 

In the course of presenting our calculations of the stopping power we make frequent 
contact with the existing literature and unify many previous results; thus this work serves 
as a review. On the other hand, we do present some original results: we introduce new 
methods, obtain an analytic expression for the stopping power in a fully ionized plasma 
that is more accurate than all previous expressions, and we provide a more precise definition 
of the Fokkcr-Planck equation for dilute plasmas. In order to make this Report accessible 
to readers who are not experts in plasma theory or have no familiarity with dimensional 
continuation or both, it is largely self-contained and is written in a pedagogical style. 
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II. METHOD 



Many physical problems involve both large and small length scales and are governed by 
a small parameter g, with < (? ^ 1. Two different physical mechanisms dominate at 
short and large distances. An infrared (IR) mechanism dominates at large distances or low 
energies; an ultraviolet (UV) mechanism dominates at short distances or high energies. In 
plasma physics, the long-distance, collective effects (described to leading order by a dielectric 
function) are the dominant infrared effects that set the long-distance scale. Hard Coulomb 
scattering, cut off by either the classical minimum approach distance or by the quantum 
maximum momentum transfer, is the ultraviolet mechanism that sets the short-distance 
scale. A novel application of dimensional continuation has recently been introduced [9] to 
treat such problems when they can be formulated in spaces of arbitrary dimensionality v. 
If the spatial dimensionality v is analytically continued below the physical z/ = 3, then the 
infrared mechanism dominates for all scales and, since it is dominant, its contribution is 
thus easy to compute in the lower spatial dimensions. On the other hand, if v is continued 
above z/ = 3, then the ultraviolet mechanism dominates for all scales, and this different 
contribution is easy to compute in the higher spatial dimensions. As a simple example, let 
us consider a case in which the dominate infrared mechanism for z/ < 3 gives the leading 
contribution of order g'^^^'^^^\ while the dominate ultraviolet mechanism for z/ > 3 gives the 
leading contribution of order ^^-(i'-s). 

^2-(-3) for 1/ > 3 (UV) , 

^2-(3-.) for z/ < 3 (IR) . (2.1) 

The actual v dependence is slightly more complicated for the problem we will study; there- 
fore, we shall look at this simple example first since it illustrates the point more concisely. 
Prom Eq. (2.1), we see that when the infrared contribution is analytically continued from 
< 3 to > 3 it becomes subleading compared to the ultraviolet contribution since 
^ gr-C'^-s) for 1/ > 3. Conversely, when continued to < 3, the ultraviolet mechanism 
becomes subleading compared to the infrared mechanism. One concludes that the sum of 
the two processes contains both the leading and the (first) subleading terms for all spatial 
dimensionality u near the physical value = 3, and hence this sum provides the correct 
result in the physical limit — > 3. 

In general, the dominant infrared mechanism will contain a pole (3 — u)"^ refiecting an 
ultraviolet divergence that appears in this mechanism when z/ — > 3 from below, while the 
dominant ultraviolet mechanism will contain a pole (i/— 3)~^ refiecting an infrared divergence 
in this mechanism when i/ — > 3 from above. Since the physical problem can be formulated 
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in arbitrary dimensionality about u — 3 with no infinities when u — 3, these two poles must 
cancel. Residues of the poles bring in logarithms of the ratio of the relevant short and long 
distance scales (or high and low energy scales), giving a stopping power of the generic form 
dE/dx — B g^lnCg^ -\- 0{g^) . Often, this ratio is large, giving a large logarithm. It must be 
emphasized that our method evaluates not only the coefficient B of such large logarithms, 
but also the constant term C underneath the logarithm. This is so because it computes both 
the leading and first subleading terms. Often in a physical problem it is easy to compute the 
large logarithm, but the computation of the constant under the logarithm cannot be done or 
is very difficult to do. The new dimensional continuation method makes this easy. Since the 
method just described is somewhat subtle, another simple but fully pedagogical example of 
how it works out is given in Appendix A. 

One could object that we do not explicitly prove that larger subleading terms are not 
present: one may ask if an additional term that has a power dependence between gr2~(''~^) 
and can appear. However, simple dimensional analysis shows that such terms of 

intermediate order cannot be formed. The point is that, in examples such as the one we 
consider, only two physical mechanisms dominate, one at large and the other at small scales. 
Since the two mechanisms involve different physics, it is natural that different combinations 
of the basic physical parameters come into play, and hence give quite different dependencies 
on the small parameter g when the dimension v departs from v — 3. 

We have illustrated the method of dimensional continuation by a very simple model of 
the i^-dependence of the coupling parameter g given in (2.1). For the case we shall examine, 
however, the dependence on v is somewhat more complex, although the same dimensional 
continuation arguments apply. For a plasma, we shall see that the leading infrared and 
ultraviolet mechanisms behave as g'^~(^~'^^ and g^, respectively. When the infrared term 
is analytically continued from i/ < 3 to > 3, it becomes subleading since g|2+(i'-3) ^ g2 
for z/ > 3. Conversely, when we analytically continue from z/ > 3 to < 3, even though 
the ultraviolet mechanism has no v dependence, it nonetheless becomes subleading since 
g^ < for u < 3. The parameter is a dimensionless constant proportional to the 

overall plasma number density n. The additional parameters needed to form a dimensionless 
coupling g'^ involve the electric charge e and temperature T. At u = 3, the coupling is of the 
generic form g^ = e^n/T^, in agreement with Eq. (2.2) below. As will be seen explicitly in 
what follows, the leading hard process for > 3 involves scattering and is thus proportional 
to n giving a g'^ contribution as stated here. The leading soft process for z/ < 3 is essentially 
j • E heating which, for dimensional reasons, is proportional to where k is the Debye 

wave number, giving a g^^~^^ — g^~^^~'^^ contribution as stated here. 
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In our case, and ai u — 3 dimensions, the small parameters are the plasma coupling 
parameters^ 

9^^ = ^^^-^,^ (2-2) 

where T5 = fi^^^ is the temperature of plasma species h measured in energy units^, Cp is the 
charge of the projectile whose energy loss is being considered, is the charge of a plasma 
species labeled by 6, and is the Debye wave number of this species, which has density n;,, 
so that 

4 = (3i,elnb . (2.3) 
The total Debye wave number of the plasma is given by 

4 = E,'^^ (2-4) 

The classical dimensionless parameter (2.2) is the ratio of the electrostatic interaction energy 
of two particles of charge Cp and a Debye length apart divided by the temperature 3], of the 
plasma species 6, with the temperature being measured in energy units (as we shall always 
do) . A parameter g of this form is the correct parameter to describe plasma effects order-by- 
order because the effects come, up to logarithmic factors, in^ integer powers of g. To make 



^In this paper we use rationalized cgs units so that, in three-dimensional space, the Coulomb 
potential energy has the form e^/(47rr). Wc do this because then no factor of 47r appears in 
Poisson's equation for the potential, a factor that is peculiar to three-dimensional space, and we 
shall need to work in a space with z/ 7^ 3 dimensions. 

^Although we shall often graph results for a plasma whose various components are at a common 
temperature, for completeness we shall work in a general case in which each plasma species h is in 
thermal equilibrium with itself at temperature T^. 

Q 2 

The proper plasma expansion parameter has the generic form 9 = (3 ^ k- Thus it is related to 
the often used plasma parameter F = /? |^ ^4|n^i/3 p3 _ g^/^_ The correct integer powers of g 
which appear in all perturbative expansions of plasma processes appear as fractional powers of T. 
That plasma perturbation theory involves integer powers of g is discussed in detail, for example, in 
Brown and Yaffe [5] in footnote 26 and in Section 3 of Appendix F. The plasma coupling parameter 
g also appears explicitly in the BBGKY equation chain if times are scaled by the inverse plasma 
frequency and lengths are scaled by the Debye length. See, for example. Section 12.5.1 of Clemmow 
and Dougherty [10] who denote g by 1/nh^. It is worth noting that the inverse of the number of 
particles in a sphere whose radius is the Debye length is given by [(47r/3)«;~^n]~-'^ = 3g. 
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an explicit (albeit slightly arbitrary) definition of the overall couphng of the projectile to the 
plasma, we define 



b 

Our calculation gives the energy loss to the generic order [In g"^ + const] = In C g^ 
in the plasma parameter, including the constant C, and to all orders in the parameters 

that measure the strength of the interaction of the projectile of charge Cp with a plasma 
particle of charge e^, with a typical or average relative velocity Vpb between the projectile 
and plasma particle. The presence of Planck's constant h in the denominator shows that 
this is a quantum-mechanical. Coulomb coupling parameter. When Vpb becomes large, fjpb 
becomes small. This is equivalent to the formal limit of large H. Hence, when fjpb is small, 
quantum effects may be important. The plasma is taken to be composed of non-relativistic 
particles that have no degeneracy so that they are described by classical, Boltzmann statis- 
tics. We show in Appendix B that the 'Convergent Kinetic Theory' method of Refs. [11], 
[12], and [8], when evaluated in the leading order in which it was derived, produces the 
same results that are produced by our method of dimensional continuation. However, that 
method generally produces spurious, higher-order corrections in the plasma parameter that 
must be discarded. They must be discarded because they do not include all of the terms 
of the given order in the plasma coupling parameter. The inclusion of these terms gives, 
in general, misleading results^. Our method has the virtue of producing only the leading 
order terms unaccompanied by any other spurious, higher-order terms, terms that must be 



striking example of how the retention of only a part of the terms in a given order can give a 
very misleading result is provided by the calculation of the energy variation of the strength of the 
strong interaction in elementary particle physics. The easiest part of the computation is to obtain 
the effects of virtual quarks, which is akin to traditional calculations in quantum electrodynamics. 
If this is done, one concludes that the strength of the interaction increases with energy. However, 
these quark terms are only a part of the leading-order result. In this same order, the contribution 
of virtual gluons overwhelms that of the quarks and the total, complete result shows that the 
interaction strength decreases as the energy increases. An even more blatant example of the error 
of principle entailed in keeping only some but not all terms of a given order is provided by the 
instruction in an elementary physics lab: the sum of 2.1 and 2.123456 is 4.2, not 4-223456; it is 
inconsistent to retain more decimal places than those of the number with least accuracy. In physics, 
half a loaf is not better than none. 
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deleted in other methods. Some authors, for example Refs. [13] and [14], retain the spurious 
higher-order terms and thus provide inconsistent results. 

Exactly the same considerations apply to the calculation of the momentum loss as a 
projectile traverses a plasma. This result, taken together with the energy loss computation, 
uniquely determine the coefficients in a Fokker-Planck equation that describes the general, 
statistical, motion of particles in the plasma. Such coefficients are sometimes described as 
"Roscnbluth potentials" which were introduced in Ref. [15] and discussed in several places, a 
good reference being Ref. [16]. Our coefficients, however, contain no arbitrary parameters or 
cutoffs and are well-defined with no ambiguity to the order g"^ In Cg^ to which we work. We 
also apply the same methods to calculate the rate of equilibration of two plasma components 
at different temperatures, again with leading and ncxt-to-leading accuracy in the plasma 
coupling g. We shall postpone the derivation and description of the Fokker-Planck equation 
until we have first presented our results for the energy loss or stopping power dE/dx for 
several cases of interest. 

We have stated that our result is generically of order g^\\.ng'^ + const]. This gives the 
correct order as far as the plasma density n is concerned, namely, discarding other parameters 
that are needed to provide the right dimensions, the result is of order n[lnn + const]. But 
we should describe the accuracy of the result in this paper with more care. It is of the form 



(2.7) 



Since ki, has the dimensions of inverse length, while e^Ki, has the dimensions of energy, the 
prefactor e^nl has the proper overall dimensions of energy per unit length. Thus the functions 
and Q are dimensionless functions of dimensionless variables — Vp is the projectile velocity 
and the velocity Vb is the average thermal velocity of the plasma species^ h. The functions 
!F and Q may also depend upon ratios of all the particle masses that are present. Since 
e^K^ — An g^jjT (T/el), the overall factor is of the generic form g'^. A point to be made is 
that this leading order calculation, which has the formal overall factor must, by simple 
dimensional analysis, involve overall factors of the dimensionless parameters gip^, which are 
of first order in the plasma density. Powers of the dimension bearing electric charge must 
appear in either a quantum-mechanical parameter 77 or in the classical plasma coupling 
parameter g, with g"^ bringing in a factor of the plasma density. Long ago, Barkas et al. [17] 



^Whcn f]pi) is small, the function Q has a term of order Inf/p^ that adds to the Ingpiy, giving 
^'^{dphMpb) = ^^iP^KiiVph) which converts the classical short-distance cutoff into a quantum cutoff. 
This is in keeping with the remark above that quantum effects may become important when 77^^ 
becomes small. 
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found differences between the ranges of positive and negative pions of the same energy. This 
imphes that there are corrections to the energy loss of cubic order in the projectile charge^, 
terms of order e^. In a plasma, such 'Barkas terms' must involve dimensionless parameters 
and are thus necessarily of order rfi/'^ in the plasma density. These Barkas terms, as well 
as other terms of order g^, are one higher order in the plasma coupling g to which we shall 
work. 

The usual method for obtaining the energy loss for a charged particle moving through 
matter is to divide the calculation into two parts: the long-distance, soft collisions and the 
short- distance, hard collisions. Collective effects are important in the long-distance part, 
and it is evaluated from the j ■ E energy loss of a particle moving in a dielectric medium. The 
hard collisions are described by Coulomb scattering. The rub is to join the disparate pieces 
together. For the case of classical scattering, this is often done by computing the energy 
loss in Coulomb scattering out to some arbitrary long-distance, maximum impact parameter 

and then adding the j • E energy loss integrated over the space outside of a cylinder 
whose radius is this maximum impact parameter. The hard scattering processes within 
the cylinder produce a logarithmic factor ln(i?/6jiuii), where 6min is the minimum classical 
distance of closest approach in the Coulomb scattering. The soft, collective effects outside 
the cylinder produce a factor involving the Debye radius k~^, ln(l//tj5i?), which has the same 
overall outside factor. Thus the arbitrary radius B cancels when the two parts are added. 
Hence such methods must yield the correct coefficient of the large logarithm Kohuan), 
and they do so without much difficulty of computation. However, the purely numerical 
constants that accompany the logarithm (which are expected to be of order one) are harder 
to compute. Here we describe an easily applied method that yields both the constants in 
front of and inside the logarithm, with no spurious higher order terms being introduced along 
the way (as is the case with other methods). The new idea is to compute the energy loss 
from Coulomb scattering over all impact parameters, but for dimensions z/ > 3 where there 
are no infrared divergences. A separate calculation of the energy loss using a generalization 
of the j ■ E heating is done for z/ < 3, where the volume integration may be extended down to 
the particle's position without encountering an ultraviolet divergence. Both of these results 
have a simple pole at z/ = 3, but they both may be analytically continued beyond their 
initial range of validity. In their respective domains, z/ > 3 and z/ < 3, both calculations 
arc performed to the leading order in the plasma density. As will be seen, although the 
Coulomb scattering result is the leading order contribution for z/ > 3, it becomes subleading 



As we shall see, corrections involving the quantum parameters r]pi, are even functions that are 
unchanged by the reflection r] — > —rj. 
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order when u < 3. Conversely, the j • E heating is subleading for i/ > 3 but leading for i/ < 3. 
Hence, the sum of the two (analytically continued) processes gives the leading and (first) 
subleading terms in the plasma density for all dimensions u, and thus, in the limit — > 3, 
the pole terms of this sum must cancel with the remainder yielding the correct physical limit 
to leading order in the plasma density. 

The highly ionized classical plasma with which we are concerned is described exactly 
by a coupled set of kinetic equations, the well-known BBGKY hierarchy as described, for 
example, in Section 3.5 of Ref. [18]. This fundamental theoretical description makes no 
explicit reference to the spatial dimensionality, and hence it is valid for a range of spatial 
dimensions u about i> — 3. Thus our method of dimensional continuation may be applied to 
a plasma. This hierarchy holds for arbitrary plasma densities. We are interested, however, 
in the computation to leading order in the plasma density of the energy loss of a particle 
traversing the plasma. The leading low-density limit of the BBGKY hierarchy changes as the 
spatial dimensionality z/ changes. For u < 3, the long-distance, collective effects dominate, 
and the equation derived by Lcnard and Balescu applies [6,7]. An alternative derivation 
of their result is presented by Duprcc [19], and a clear pedagogical discussion appears in 
Nicholson [20]. Clemmow and Dougherty [10] provide a derivation of the Lenard-Balescu 
equation and prove that it shares the basic features of the Boltzmann equation; namely that 
it conserves particle number, total momentum and energy, and that it obeys an H-theorem 
(entropy increases) like the Boltzmann equation with the long-time, equilibrium solution 
being a Maxwell-Boltzmann distribution. The equation describes the interaction of the 
various species that the plasma may contain. In the limit in which one species is very dilute, 
as is our case in which we examine the motion of a single "test particle" or "projectile" moving 
through the plasma, the energy lost in the particle motion is described by a generalization 
of its j • E Joule heating with the background plasma response given by the permittivity of 
a collisionless plasma. On the other hand, when the spatial dimension u is greater than 3, 
the short-distance, hard Coulomb collisions dominate. For these dimensions, the leading low 
density limit of the BBGKY hierarchy is described by the familiar Boltzmann equation.^ The 
Boltzmann equation is derived, for example, in Section 16 of Ref. [21] and also in Section 3.5 
of Ref. [18]. We use the Boltzmann equation to obtain the leading order energy loss rate 
when u > 3. Since we are concerned with the motion of a single "projectile", the Boltzmann 
equation reduces to the product of the energy loss weighted cross section times the plasma 



In this case, we may go beyond the classical BBGKY hierarchy limit in that we may use the full 
quantum-mechanical cross section rather than its classical limit in the Boltzmann equation. The 
validity of this extension is, however, obvious on physical grounds. 
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density. The derivations that we have just described, which start from first principles, justify 
the methods outhned in the previous paragraph, the methods that we shall use. 

In Ref. [9] , the method was illustrated by the simplified case in which the charged particle 
moved through a dilute plasma with a speed that is much larger than the speeds of the 
thermal electrons in the plasma. Here we shall extend that work to the case in which the 
charged particle projectile moves with arbitrary speeds. As we have noted above, we work 
in the dilute limit in which the plasma density is our small parameter.^ It is obvious that 
the energy loss for i/ > 3 as calculated for scattering is proportional to the first power of 
the plasma density, dE^/dx ~ n. On the other hand, as we shall see explicitly below, the 
computation of the energy loss for i/ < 3 behaves as dE^ /dx ~ Thus we see 

explicitly that the infrared computation that accounts for the collective effects in the plasma 
and gives the leading term for i/ < 3 becomes non-leading when it is analytically continued to 
V > 2). Conversely, the ultraviolet, hard scattering computation, which is leading for i/ > 3, 
becomes subleading when analytically continued to v < We conclude that the leading 
order (in density) energy loss of a projectile particle moving in three dimensions is given by 

dE {dE^- dE<} 

= lim I — — + — — } . (2.8) 



dx i^-*3 \^ dx dx ^ 

To the order in the coupling g to which wc work, a plasma species b is described completely 
by a phase space density /^(r, p, t). The projectile p may also be described by a phase space 
density fp{r, p, t) that contains delta functions restricting the momenta p and the coordinates 
r to be those of the projectile's trajectory. The projectile energy loss (or "stopping power") 
that we deal with is the rate of kinetic energy loss, the time derivative of (in u dimensions) 

E^^ f d'r f -^^^ f (r,p,t). (2.9) 

The total kinetic energy of plasma particles of species b is given by 

The Lenard-Balescu and Boltzmann equations that we use to derive the projectiles' energy 
loss obey an energy conservation law that entails only these kinetic energies. 



^One should work with a dimensionless parameter so that stating that it is small is unambigu- 
ous. For our case, the dimensionless parameter is the square of the plasma coupling parameter, 
^2 = (e^Ko/'iTrTf = (eVlSTT^T^) n, but to save writing we shall simply use the density n as our 
parameter. 
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There are, of course, additional contributions to the total energy — potential energy con- 
tributions that involve collective plasma effects. The kinetic energies are of zeroth order in 
the plasma coupling g and their time derivatives, which come about because of the Coulomb 
forces, are of order g'^ In g'^. On the other hand, the potential energies (and possible collec- 
tive effects) are of higher order in the coupling g and their time derivatives are of an order 
that is higher than that to which we compute. Hence these potential energy terms do not 
contribute to the energy balances accounted for by the Lenard-Balescu and Boltzmann equa- 
tions, and only the kinetic energies are relevant. We conclude that, to the accuracy to which 
we compute, we have an unambiguous partition of the projectile's energy loss into energies 
gained by individual particle species in the plasma. We define the energy loss dE/dx of the 
projectile to be positive, 



dE _ 1 dEp 
dx Vp dt 



(2.12) 



and we have 

dx dx 



y ^ , (2.13) 



where 

dEb _ 1 dEi, 
dx v„ dt 



(2.14) 



defines the energy loss of the projectile p to the plasma particles of species b or, equivalently, 
the energy gain of the plasma particles b brought about by the projectile p moving through 
the plasma with velocity Vp. Such an unambiguous partition into the energy gained by the 
individual species in the plasma does not hold in higher orders, and so such an accounting 
cannot be done for strongly coupled plasmas which are entangled with collective excitations. 

An important check of the validity of our results is that they satisfy the condition that 
the total energy loss vanishes for a swarm of projectile particles in thermal equilibrium with 
the plasma through which it moves. We must have 

dE\ (I3m,\"^ /- 3 r mpvl\ dE 

'rfvpexp<^ — -—\vp—^Q. (2.15) 



dt / \2n J J f ^ [ 2T } ^ dx 

Thus the energy loss dE/dx must become negative for low projectile velocities Vp so that 
its integral over all velocities can vanish. As we shall see, this constraint is an automatic 
consequence of the method that we employ. 
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We turn now to describe the results of the calculation of dE/dx. In addition to parti- 
tioning the energy loss into different particle species that make up the plasma, we shall also 
increase the generality of our results by assuming that, although each species h is internally 
in thermal equilibrium, it may have a private temperature which differs from species to 
species. We shall briefly compare our results in the classical limit with PIC simulation data^, 
illustrate our stopping power under solar conditions, and provide a more lengthy exposition 
relevant for inertially confined laser fusion experiments. These examples illustrate the cal- 
culation of dE/dx within our method, and we can directly compare our results for dE/dx 
with those that are typical of the current literature. The dimensional continuation method 
that we have developed and applied to the calculation of the stopping power can be used 
to calculate other physical processes. As a final example, we also will use this method to 
compute the rate at which Coulomb interactions in a dilute plasma bring two species into 
thermal equilibrium. The mass ratios and initial temperatures are arbitrary, and quantum 
corrections are included. Like the stopping power, the temperature equilibration rate^° is 
of the form F = Bg^lnCg'^ -\- 0{g^), and we calculate the prefactor B and the constant C 
under the logarithm exactly. After presenting these results, we shall describe the general 
formulation that results in a Fokker-Planck equation. The range of vahdity of the Fokker- 
Planck equation will be assessed, and then we shall present the derivations of details of our 
results. 

III. RESULTS 

For z/ > 3, the BBGKY hierarchy reduces, in our dilute plasma limit, to the Boltz- 
mann equation with classical Coulomb scattering. In this higher-dimensional space, only the 
squared momentum transfer weighted Coulomb collision cross section enters into the rate 
of energy loss. Although quantum-mechanical corrections arc not significant in the lower- 
dimensional z/ < 3 spatial regions where long-distance effects dominate, they may be impor- 
tant in the higher-dimensional > 3 space in which the short-distance collisions dominate. 



^In the particle-in-cell (PIC) technique as applied to plasmas, the plasma species are modeled as 
collections of quasi-particles, each representing a large number of real particles, moving through 
a numerical grid. The particle positions and a weighting factor are used to assign electric charge 
densities to the nodes of the grid. Poisson's equation is then solved for the potentials at the 
nodes, which gives the electric fields at the nodes, and finally an inverse weighting factor is used to 
determine the electric field at the particle positions. The particles are then moved using Newton's 
equations to start the next time step. 

^'^It should be clear from the context in which it is used, whether the letter F stands for a rate 
rather than the traditional plasma coupling mentioned in footnote 3. 
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And although our approach is initially based on the classical BBGKY hierarchy, it is physi- 
cally obvious that quantum corrections must be incorporated in the higher-dimensional re- 
gion when they become important. The dimensionless parameter that distinguishes whether 
or not quantum effects must be taken into account is (in the i/ — > 3 limit) 



Vpb 



(3.1) 



ATrhvpfj ' 

where Vpb is a typical relative velocity of the projectile (p) of charge and a plasma particle 
(6) of charge e^- The limit of large fjpb describes slow particles. This limit is equivalent to 
the formal /i — > limit, and thus the classical calculation applies here. However, when fjpb is 
not large, quantum effects must be taken into account for v > 3. We shall first examine the 
classical case and then later the quantum corrections to it. 



A. Classical Regime 

The classical results to the order in the plasma coupling to which we compute are sum- 
marized in Sec. IX. The complete energy loss to the plasma species b in the classical case is 
given by 



(3.2) 



dx dx dx 

where the two contributions are contained in Eq's. (9.5) and (9.7), and with the aid of 
Eq. (4.21) they can be written as:^-*^ 

.2 ..2 / \ 1/2 



dx 



4:71 rUpVp 



rrib 



In A 



27r/?5 

CpCb K nib u 



J duu^''^ exp | — -/Sfem^fp M j> 



47r rupb 1 — u 



+ 2-27 



PbMpbv;-- 

^ u 



(3.3) 



where 7 ~ 0.5772 is Euler's constant, and 



dEt 



dx 



IE. 




1: 


4:71 
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1 



Ptotal(^^pCOS( 



-F{vpCOS 9) In 



471 27r PbUlpV^ ptota\{Vp) 



F*{vp) In 



F{vp cos ( 



(3.4) 



^^To save writing, we use e to denote the absolute value of the charge of a particle. Thus CpCb is 
always positive even if projectile (p) and plasma (6) particles have charges of opposite sign. 
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Here K is an arbitrary wave number. As we shall soon show, the total result (3.2) does 
not depend upon K. However, sometimes choosing X to be a suitable multiple of the Debye 
wave number of the plasma simplifies the formula. 

We use Vp to denote the speed of the projectile of charge Cp and mass rup whose rate of 
energy loss in the plasma we are computing. Rationalized units are used for the charge so 
that, for example, the Coulomb potential energy in three dimensions reads ep/(47rr). We 
write the inverse temperature of the plasma species b as Pi, — T^^, which we measure in 
energy units. The charge and mass of the plasma particle of species h are written as e;, and 
mb, with the corresponding Debye wave number k^, of this species defined by 

nl = P^el rib , (3.5) 

where rib is the number density of species b. The total Debye wave number is defined by 
the sum over all the species 



^ = E.'^^ (3-6) 



The relative mass of the projectile and plasma particles is denoted by rUpb, with 

' (3.7) 



rupb rrip rUb 
while 

Mpb = rrip + mb (3.8) 

is the corresponding total mass. 

The function F{u) is related to the leading-order plasma dielectric susceptibility. As in 
the discussion of Eq. (7.8), it may be expressed in the dispersion form 

F{u) = - dv ^^"^"'^""^ , (3.9) 
i_oo u-v + iri 

where the limit 77 ^ 0"*" is understood. The spectral weight is defined by 

Ptotal(^^) = Pc{v) , (3.10) 

where 



Pc(^^) ^/^c^Y^^expl-^/^cmct;^! . (3.11) 

It is worthwhile noting here several properties of F that will be needed throughout for an 
understanding of the results. Clearly the spectral weight is an odd function. 



16 



Pc{-V) = -Pc{v) . 



(3.12) 



Hence the variable change v ^ —v in the dispersion relation (3.9) gives the reflection property 

F{-u)=F*{u). (3.13) 

These properties imply that the total integral in the lower-dimensional contribution (3.4) is 
real as it must be. Since 

Im^— ^7rS(x) , (3.14) 

X — IT) 

the dispersion form (3.9) gives 

F{u)-F*{u)^27riptot.i{u) 

= F(u) - F{-u) , (3.15) 

with the second equality just a repetition of the reflection property (3.13). These results 
show that 

To show that dE/dx is independent of the wavcnumbcr wc take the analogous deriva- 
tive of Eq.(3.3) and make the variable change u = cos^O. One term involves duu^/'^ = 
2d cos 9 cos'^9, the other duu^^'^u^^ = 2d cos 9. For the integral involving this second term, 
we insert 1 = dcos9/dcos9 in the integrand, and integrate the result by parts. In this way 
we find that the total result (3.2) is indeed independent of the arbitrary wave number K. 

The classical result applies for a low velocity projectile moving in a relatively cool plasma. 
In addition to such physical applications, our results, which are rigorous and model inde- 
pendent to the order g"^ In Cg^ to which we compute, may be used to check the validity of 
computer calculations such as those utilizing classical molecular dynamics, as illustrated in 
Figs. 1 and 2. Such calculations must agree with our results in those regions where the 
plasma coupling g is not large. The following figures display our analytic results for the 
energy loss (3.2) using the classical expressions (3.3) and (3.4). They arc compared with 
the molecular dynamics calculations of Zwicknagel, Toeppfer, and Reinhard [22], as cited in 
Ref. [14]. See also the review of Zwicknagel, Toeppfer, and Reinhard 
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FIG. 1. Comparison of our classical result dE'^/dx (curve) versus projectile velocity to PIC simulation 
data of Zwicknagel et al. cited in Ref. [14] (squares). The projectile is a very massive ion of charge Z = b 
moving through an electron plasma with n = 1.1 X lO^'^cm"'^ at T = 14 eV = 1.6 X 10^ K, giving a plasma 
coupling gp = 0.21. The thermal velocity of the electron Vf, = y'ST/m^ = 2.6 X 10^ cm/s has been used 
to set the velocity scale. 
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FIG. 2. Comparison of our classical result dE'^/dx (curve) versus projectile velocity to PIC simulation 
data of Zwicknagel et al. cited in Ref. [14] (squares) for a very massive ion of charge Z = 10 moving through 
an electron plasma with n = 1.4 X 10^'^ cm"'' at T = 11 eV = 1.3 X 10^ K, giving a plasma coupling 
Qp = 0.61. The thermal velocity of the electron is = 2.4 X 10^ cm/s. 

B. Quantum Corrections 

Thus far, our discussion has been for those cases in which classical physics apphes. In 
these cases, the quantum parameters defined in (3.1), 
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(3.17) 



4:TTh Vpb 

are large. In the energy loss problem, these are the only independent dimensionless param- 
eters that entail the quantum unit, Planck's constant H. The parameters are large when the 
average relative velocity Vpb is small which, as far as an fjph parameter is concerned, corre- 
sponds to the formal hmit — > 0. In this section, we treat the general case where the size 
of the quantum parameters fjpf, has no restriction. 

According to our dimensional continuation method and, in particular, the discussion 
leading to (10.27), the general case is obtained by adding a correction to the classical result 
(9.6). Namely, the energy loss to the plasma species b in the general case appears as 

dEb^dEl^dEl^ 
dx dx dx ' 

where, we recall, the classical contribution dE^/dx is described by Eq.'s (3.2) - (3.4), while 
Eq's. (10.27) and (4.21) give 



dE^ 




/ (3bmb 


dx 


47r 2[5bmpvl 


V 2n 
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Mpb Vp 






nib Vpb 






1 + 


Mpb Vp 






rrib Vpb 



1/2 



dv, 



'pb 



cxp 



+ 1 



exp <j -^pbirib {vp + Vpbf 



{2ReV'(l + i77p6) -lnr?JJ 



(3.19) 



Here 



npb = 



CpCb 



(3.20) 
(3.21) 



4:TThVpb 

and ip{z) is the logarithmic derivative of the gamma function. As explained in the derivation 
of (10.17), we may write 



Reip{l + if]) 



77^ 



i^kk"^ -\- rf 



-7 



(3.22) 



To illustrate our general results, we now present several plots of the total energy loss 
dE/dx as a function of the incident projectile velocity Vp. In order to have some comparison 
with previous work, we have chosen to plot the results of Li and Petrasso [24] together with 
our results. 
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Li and Petrasso [24] have evaluated the energy loss dE/dx to leading order in the plasma 
coupling, the order In contribution, and they have added an estimation of the constant 

under the logarithm, the term of order g"^ InC that we compute exactly. They do this by 
working with a Fokker-Planck approximation to the Boltzmann equation. They then define 
a Coulomb logarithm by using a somewhat arbitrary minimal classical impact parameter 
that is then corrected in a rather ad-hoc fashion to take account of quantum-mechanical 
corrections. A term involving a step function is added to the formula to correct for long- 
distance collective effects. Using our notation, the final result of Li and Petrasso [24] appears 



as 



dE,, 



e; 1 ^ 



dx Air Vp PiiTJib 



Here 



where 



and 



^ mb d 
nip dy 



H{y) = ^ 

nip 



'71 



1 

dy^ 



dzz''^ 



M + 9{y - 1) In (2e-^yV2) ^ 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



with 6{x) the unit step function: 6{x) = for x < and 6{x) = 1 for x > 0. Li and Petrasso 
[24] define a Coulomb logarithm in terms of the combination of classical and quantum cutoffs 
as described above, namely 



In Ah = — - In Kr. B, 



where 



\2nipbUbJ \A7rnipbulJ 



+ 



(. 



b > 



€p6b 



Y 



(3.27) 



(3.28) 



(3.29) 



in which nipb is the reduced mass of the projectile (p) - plasma particle (6) system, and 

2 

defines an average of the squared projectile and thermal velocities. 

Figures 3-5 compare the results of Li and Petrasso [24] with our results for a proton 
projectile in a fully ionized hydrogen plasma, a neutral plasma of electrons and protons, 
with the electrons and protons at a common temperature — Tp — T. These comparisons 
are made over plasma temperatures and densities where the plasma coupling parameter g is 
reasonably small so that our approximation is essentially exact. 
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FIG. 3. Proton projectile in a fully ionized hydrogen plasma: dE/dx (in MeV//ini) us. Vp/vg. The solid 
line is the result from this work (BPS), and the dashed line is the result of Li and Petrasso (LP). Here 
the temperature T = 1 kcV and density rie = 5 X 10^^ cm"'' are chosen to correspond approximately to 
values at the core of the sun. The plasma coupling is Qp = 0.061 and the thermal speed of the electron is 
Vg = ^J^TjuLe = 2.30 X 10^ cm/s. The electron fugacity Ze = exp{/3//e} has the value Ze = 0.25, and so 
the relative Fermi-Dirac statistics corrections are Zf.jl^l'^ ^ 9% for both LP and BPS. 
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FIG. 4. Proton projectile in a fully ionized hydrogen plasma: dE/dx (in MeV//xm) vs. Vp/ve with 
T = 0.2 keV, He = 10^^ cm"'', giving a plasma coupling gp = 0.097, and the thermal speed of the electron 
Ve = 1.03 X 10*^ cm/s. The solid line is the result from this work (BPS), while the dashed line is the result 
of Li and Petrasso (LP). 
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FIG. 5. Same as in Fig. 4 except that the temperature is raised to T = lOkeV, giving Qp = 0.00027 
and Ve = 7.26 x 10^ cm/s. There is less of a discrepancy between the results of Li and Petrasso (LP) and 
our work (BPS) at higher temperatures. 



The figures show fair agreement between our results and those of Li and Petrasso in regions 
where the projectile velocities are large, but significant discrepancies appear for small and 
intermediate projectile velocities. These discrepancies are related to the fact that the rate of 
energy loss for Li and Petrasso (LP) does not keep a swarm of projectiles in thermal equilib- 
rium with the plasma particles even though the swarm has the same initial temperature as 
that in the plasma, whereas we will show in Eq. (4.26) that our (BPS) energy loss expression 
does maintain thermal equilibrium. This is to say, the thermal average for Li and Petrasso, 
given by 

3/2 

dx 



dt \27: ) J dx 



(3.30) 



does not vanish, while the thermal average of our energy loss does vanish 

dEm 



dt 



, (3.31) 



as we will show in Eq. (4.26) below. 

It is easy to prove that (3.30) does not vanish. We write 



dE, p \ / dE, p \ / dEj „ , . . 

+ ( —7^ ) , (3.32) 



dt / ^ \ dt / ^ \ dt 



G 



corresponding to the separate contributions involving the functions G{y) and H{y) in formula 
(3.23) of Li and Petrasso. Clearly 
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1 + 



(3.33) 



Since 2e — 1.1229 • • •, the additional logarithmic contribution in (3.26) (with y > 1 because 
of the step function factor) is also non- negative. Hence H{y) > 0, and we conclude that 



dt 



> 0. 



(3.34) 



The form (3.24) of G{y) gives 



d 



exp <i --pnipvl \ /i ( ^/Snibvl 



(3.35) 



(5mpVp dvp 

Thus evaluating the velocity integral in spherical coordinates and integrating by parts yields 



dE^ 
dt 
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(3.36) 



Here 



— InAb 
dv„ 



1 dBl 



2Bl dVp 



>0, 



2mpbUbJ \mpbul) 



ut 



and hence we have^^ 



dE^ 
dt 



> 0. 



(3.37) 



(3.38) 



Thus the formula of Li and Petrasso violates the thermal equilibrium condition, a condition 
that our results always obey. 



^^Since fi{y) vanishes when y ^ and is finite when y — ^ oo, there are no end-point contributions 
in the partial integration and the resulting integral is well defined. 
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FIG. 6. Proton projectile moving in a fully ionized hydrogen plasma. The classical and quantum contri- 
butions to the BPS stopping power dE/dx of this work (in MeV//J.m) are shown vs. Vp/ve- As in Fig. 3, 
T = 1 keV and rig = 5 x 

1025 

cm ^, giving Qp = 0.061 and Ve = 2.30 x 10^ cm/s. The solid hne is the 
purely classical result (3.2) summed over species, dE'^ /dx, while the dashed line is the quantum correction 
(3.19) summed over species, dE'^/dx. 

C. Quantum vs. Classical Contributions 

Figure 6 illustrates the size of the quantum corrections given in Eq. (3.19) relative to 
the classical formula (3.2). For the parameters listed in its caption, the quantum correction 
is about 40% of the classical contribution for Vp/ve greater than 2. For lesser values of the 
projectile velocity, the relative importance of the quantum correction decreases. 

D. Quantum Limit 

We started our discussion of the stopping power by examining the low velocity limit. 
This limit is contained in the classical result (3.2), which adds the regular, long-distance 
contribution (3.4) to the well-behaved sum of singular short- and long-distance contributions 
(3.3). The long-distance, collective effects are always described by the classical dielectric 
properties of the plasma. The short-distance effects presented in Eq. (3.3) involve a classical 
description of the Coulomb scattering which gives a classical minimum approach distance 
bo = {^p^b/4:n){l/mpbVpi^) that provides the short-distance cutoff. This classical description 
of the scattering is valid in those situations where the projectile velocity Vp is sufficiently 
low. The momentum transfer defines a quantum wave length of order h/{mpf,Vpf,). When 
this length is on the order of, or larger than, the classical minimum approach distance bo, 
then the detailed quantum-mechanical treatment that we have given is needed. When the 
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quantum wave length is much greater than the classical minimum approach distance 60, then 
the full quantum description of the short-distance scattering simplifies, with the first Born 
approximation sufficing. This happens when h/{m.pbVpb) ^ {epeb/4:7r){l/mpbVpfj), or when 
Vpb — {epeb/ATrhvpb) <^ 1. The stopping power in this extreme quantum limit is given by 



77 < 1 : 
dEb _ dEb^Q 



dx 



dx 



+ 



dEb,n 
dx 



(3.39) 



The regular, long-distance contribution dE^^/dx is the classical result (3.4) that is not 
changed. The dEb,Q/dx contribution is obtained by adding the large velocity fimit (10.38) 
of the difference between the total and the classical hard scattering contributions (1 e. the 
purely quantum correction) to the classical result (9.5). This gives 
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(3.40) 



This result holds when {epCb / 4:n) / (hvp) is much less than one, but there is no restriction 
on the comparison of the energy of the projectile with the temperature. When rribVp ^ T5 for 
all the plasma species b, the results simplify considerably. The limit is characterized as the 
formal Vp ^ 00 limit. In this case, the contribution of the electrons in the plasma dominate 
and the limit (7.41) found below gives 



-> 00 
dE< 



dE< 



In 



K'^f3emevl 



dx dx An 2(5e.meVp 
Adding this to the corresponding limit of Eq. (3.40) yields^^ 



(3.41) 



^^The factor fibrribVp in the exponent of Eq. (3.40) restricts the contribution of the u-integration 
to values less than or of the order of {f5bmbVp)^^ which, for large Vp, is much less than one. Hence 
the factor 1 — in the logarithm may be replaced by unity and, moreover, the contribution of 
the small-particle-mass electronic component of the plasma dominates. Since the exponent damps 
out the contributions of large u, the upper limit of the integral may be extended from u = 1 to 
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This formula is simplified and its nature clarified if we introduce the electron plasma fre- 
quency uje defined by 



^l-"^-^. (3.43) 
me Perrie 



for we now have 
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This well-known high-velocity limit is valid when the projectile velocity is much larger 
than the thermal velocity of the electrons in the plasma, Vp ^ Ve and, in addition, when 
the projectile velocity is sufficiently large that the quantum Coulomb parameter is small, 
Vp = {epCe/AnhVp) -C 1. Not that this high- velocity result is independent of the temperature 
of the plasma. 

E. Results Relevant for Laser Fusion 

We turn now to examine cases that are relevant to the deuterium-tritium (DT) plasmas 
in laser fusion capsules. In an inertial confinement fusion (ICF) capsule filled with DT gas, 
an a particle of energy Eq — 3.54 MeV is created at threshold in the reaction D + T ^ a + n. 
This a particle slows down and eventually deposits its energy into the plasma, if the range 
is short enough, or exits the ICF capsule entirely, if the range is too long. The more energy 
deposited into the plasma by the a particle then the hotter the plasma becomes, and this 



u —>■ oo. Writing the last term in the second square brackets in Eq. (3.40) as u"^/^ = 2{d/du) v}/'^ 
and integrating by parts produces no end-point contributions and cancels many terms save for one 
involving PerUpVp. Finally, changing variables to z = l3e'meVpU/2 with the aid of the integrals 

POO 

/ dzz^/'^e-' = r{3/2) = ^/^/2, 
Jo 

and 

oo 

dz In z = r(3/2) V(3/2) = r(3/2) [2 - In 4 - 7] , 
gives the limit that leads to the result (3.42). 
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in turn increases the rate of DT fusion. Obviously then, the precise value of the a particle 
range can have a dramatic impact on ICF performance. 

Again we shall assume that all the electrons and ions are at a common temperature T. 
As we shall see, our results can differ by 20% or so from those of Li and Petrasso [24], results 
that have been used in the description of such laser fusion experiments^^. We first plot in 
Fig. 7 the total energy loss to all species dE/dx for an alpha particle moving through a DT 
plasma with a temperature T = 3 keV and electron density ng = 10^^ cm~^. We shall also 
exhibit in Fig. 8 the separate energy losses to the electrons and to the ions that are composed 
of equal numbers of deuterons and tritons. Again our results are compared to those of Li 
and Petrasso. 
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FIG. 7. Alpha particle projectile traversing an equal molal DT plasma. The stopping power dE/dx 
(in MeV//im) is plotted vs. energy E (in Mev). The solid hne is the result from this work (BPS), 
while the dashed line is the result of Li and Petrasso (LP). The energy domain lies between zero and 
the a particle energy Eq = 3.54 MeV produced in the DT reaction. The plasma temperature is 
T = 3 keV, the electron number density is rig = 1.0 X 10^^ cm"'', with deuterium-tritium number densities 
f^d = fit = 0.5 X 10^^ cm"'' (for charge neutrality). The plasma coupling is gp = 0.011, and the thermal 
speed of the electron is Ve = 3.98 X lO'^ cm/s. The BPS result is essentially exact since the plasma coupling 
is so small. 



^'^N. M. Hoffman and C. L. Lee [25] have used the stopping power computations of Li and Petrasso 
to model the implosion of laser driven fusion capsules 
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FIG. 8. The results of Fig. 7 split into separate ion (peaked curves) and electron components (softly 
increasing curves). 



As shown in Fig. 8, for ion projectiles, the energy loss to the electrons in the plasma 
dominates over that to the ions when the projectile energy becomes sufficiently large on 
the scale of the temperature T. Here we provide an estimate of the energy at which this 
cross over takes place, an estimate that is valid to logarithmic accuracy, which holds when 
the logarithmic term in the energy loss formula - the "Coulomb logarithm" - is large and 
dominant. We do this for a plasma whose various species are all at the same temperature 
T. Denoting this logarithmic term by L, which we treat as a constant since its variation 
within an integral is small, the logarithmic contribution to the stopping power is contained 
in Eq. (9.5) and reads 
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(3.45) 



As we shall find, near the cross over region the projectile energy E = mpVp/2 is large in 
comparison with the temperature T, and the factor /SnifyVp in the exponent in Eq. (3.45) is 
large for ions of mass nih ~ nip. Hence only the small u region of the integration makes a 
significant contribution, and the integration region < m < 1 can be extended to < m < oo 
to obtain the leading piece. The resulting Gaussian integrals are readily done, and one finds 
that for an ion b, 
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(3.46) 
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Here we have used the ionic plasma frequency, cul — ef rib/mb and — (5mi,ujl, in the second 
equality to emphasize that the result is independent of the temperature. The cross over point 
is at a projectile energy such that PnieVp <^ 1, as we shall soon find. Hence for the electrons 
in the plasma, when 6 = e, we may approximate the exponential by unity in Eq. (3.45) to 
obtain 
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(3.47) 



On comparing these two equations, we find that the electron contribution dominates over 
the ionic contribution of species 6, 



di^c dec 
when the projectile energy E > El, with 



(3.48) 



Er 



E 



16 ruen^ Vtl^ "^6 

^. byte 



1/3 



T. 



(3.49) 



Note that the expression in square brackets here is independent of the temperature so that 
El scales hnearly with the temperature. For the parameters of Fig. 8, the estimate provided 
by Eq. (3.49) gives the cross over energy El — 0.10 MeV compared to the actual cross over 
energy of 0.18 MeV. 

The amount of energy Ei that the slowing particle with an initial energy Eq transfers to 
the ions may be expressed as 
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where 
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all ions b 

The corresponding energy loss to the electrons in the plasma is, of course, just E^ — Eq — E,. 
We can use the rough logarithmic approximations of the previous paragraph to estimate 
the energy transfer. We use the sum of the approximate ionic stopping powers (3.46) in 
the numerator of Eq. (3.50) and add the approximate electronic part (3.47) to this for the 
denominator in Eq. (3.50). Changing the integration variable E to an appropriately scaled 
velocity then yields, in this logarithmic approximation. 



EiL — 2El 



f 

Jo 



xdx 



x^ 



(3.52) 
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where El is the estimate (3.49) of the cross over energy. Since this integral damps out at 
large x values, the simple upper bound to this approximation for the ionic energy transfer 
obtained by extending the upper limit of the integral to infinity should not be too far off. A 
glance at Fig. 8 shows that the ionic energy loss is indeed dominated by small energies. We 
use 

xdx 2n 



to obtain 



1 + x^ 



3^3 



3v^' 



(3.53) 



(3.54) 



For the parameters of Fig. 8, this crude limit gives 0.24 MeV to be compared with the value 
of El = 0.38 MeV that comes from a numerical evaluation of Eq. (3.50) using our complete 
energy loss formulas with an initial energy of Eq = 3.54 MeV [c.f. Fig. 10]. In this case, the 
energy loss to the ions is small in comparison to that lost to the electrons, E^ = 3.54 — 0.38 = 
3.16 Mev. It is interesting to note that, since the approximate limit scales linearly with the 
temperature, if the plasma temperature is increased by an order of magnitude, from 3 keV 
to 30 keV, the limit of 0.24 MeV moves to 2.4 Mev. A numerical evaluation similar to that 
reported in Fig. 10 for this increased temperature of 30 keV gives E^ = 1.8 MeV. This is 
now comparable to the energy transfer to the electrons. Eg = 3.5 — 1.8 = 1.7 MeV. We 
should note that, as the temperature is increased with a fixed initial projectile energy Eq, 
the upper integration limit ^/EqJEl in Eq. (3.52) is reduced, and so its replacement by the 
limit X = oo gives an increasingly worse result. 
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FIG. 9. The distance x{E;Eq) computed by Eq. (3.55) defines an energy E{x) for a particle tliat has 
traveled a distance x starting at x = with an initial energy E{0) = Eq. An a particle of energy 
£^0 = 3.54 MeV is created from threshold in the reaction D + T — > a + n. For the DT plasma defined in 
Fig. 7, we have plotted the a particle energy (in MeV) vs. the distance traveled (in /xm). The solid line is 
the result from this work (BPS), while the dashed line is the result of Li and Petrasso (LP). They give the 
respective ranges -Rbps = 30 fim and i?Lp = 25 jim , about a 20% difference. 
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From the results shown in Fig. 7, we can compute the distance x that a projectile, starting 
with energy Eq, travels to be slowed down to reach the energy E: 

fEo f dE\~^ 
xiE;Eo) = dE i—j . (3.55) 

Figure 9 shows the inverse function, E vs. x{E;Eq), for an alpha particle with an initial 
energy Eq = 3.54 MeV, corresponding to the alpha particle produced in DT fusion. In Fig. 8 
we illustrated the energy dependence of the electron and ion components of the stopping 
power. In Fig. 10 we plot the electron and ion components as a function of the the distance 
X that the a particle has traversed. 




(3.56) 



10 20 30 

X (yU.m) 



FIG. 10. The a particle dE{x)/dx (in Mev//um) vs. x (in /im) split into separate ion (peaked curves) 
and electron components (softly decreasing curves). The energy used to compute dE{x)/dx is determined 
from the results shown in Fig. 9 while the corresponding dE{x)/dx is given by the results in Fig. 8. Again, 
the solid line is the result from this work (BPS), and the dashed line is the result of Li and Petrasso (LP). 
The area under each curve gives the corresponding energy partition into electrons and ions for this work and 
that of Li and Petrasso. For our results, the total energy deposited into electrons is E^^^ = 3.16 MeV and 
into ions is Ef^^ = 0.38 MeV, while LP gives E^^ = 3.09 MeV and El^^ = 0.45 MeV. Note that these 
energies sum to the initial a particle energy of Eq = 3.54 MeV. 

It is worthwhile comparing the results that we have just illustrated for DT produced 
alpha particles moving in a plasma at 3.0 keV and electron number density 10^^ cm~'^ with 
DT alphas moving in a hotter, more dense plasma, a plasma at 30 keV and an electron 
number density lO^^cm"^; see Figs. 11-14. In the previous case, most of the alpha particle 
energy was deposited into electrons. In the new case, much of this energy is now transfered 
directly into the ions. The new, much denser case clearly has a much shorter alpha particle 
range. 
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FIG. 11. The stopping power dE/dx (in MeV//im) of an a particle projectile traversing an equal molal 
DT plasma as a function of the projectile energy E (in MeV). The energy domain lies between zero and the a 
particle energy Eq = 3.54 MeV produced in DT reaction. The plasma temperature is T = 30keV and the 
electron number density is Uf, = 1.0 X 10^'^ cm""^, which is characteristic of plasmas for inertial confinement 
fusion shortly after ignition. The deuterium-tritium number densities are Ud = nt = 0.5 X 10^'^ cm"'' 
(for charge neutrality). The solid line is the result from this work (BPS), and the dashed line is the result 
of Li and Petrasso (LP). The plasma coupling is gp = 0.0033, and the thermal speed of the electron is 
Ve = 1-26 X 10^*^ cm/s. The BPS result is essentially exact since the plasma coupling is so small. 
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FIG. 12. The results of Fig. 11 split into separate ion and electron contributions. 
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FIG. 13. For the DT plasma defined in Fig. 11, the energy (in MeV) as a function of the distance traveled 
(in /xm) is shown for an a particle created at Eq = 3.54 MeV in the reaction D + T ^ a + n. The solid 
line is the result from this work (BPS), and the dashed line is the result of Li and Petrasso (LP). They give 
the respective ranges i?BPs = 3.7 /im and i?Lp = 2.9 /im, almost a 30% difference. 
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FIG. 14. Alpha particle dE{x)/dx (in MeV/^m) vs. x (in iim) split into separate ion and electron 
components. The energy used to compute dE{x)/dx is determined from the results shown in Fig. 13 while 
the corresponding dE{x)/dx is given by the results in Fig. 12. The solid line is the result from this work 
(BPS), and the dashed line is the result of Li and Petrasso (LP). For our result, the energy deposited into 
electrons is Ef^ = 1.51 MeV and into ions is Ef^^ = 2.00 MeV, while LP gives E^'' = 1.36 MeV and 
E^^ = 2.15 MeV. Note that both LP and BPS sum to 3.51 MeV, which is within 1% of the initial a particle 
energy Eq = 3.54 MeV. 



We conclude this section by plotting similar figures for a triton moving through a deu- 
terium plasma with T = 0.5 keV and an electron density Ue = lO^^cm"^; see Figs. 15-17. 
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FIG. 15. For a triton projectile traversing a deuterium plasma, the stopping power dE/dx (in MeV//im) 
is plotted vs. energy E (in MeV). The solid line is the result from this work (BPS), and the the dashed 
line is the result of Li and Petrasso (LP). The energy domain lies between zero and the triton energy 
Eq = 1.01 MeV produced in the DD reaction. The plasma temperature is T = 0.5 keV, the electron 
number density is rig = 1.0 x 10 cm~ , with a deuterium number density Ud = 1.0 X 10 cm~ 
(for charge neutrality). The plasma coupling is gp = 0.025, and the thermal speed of the electron is 
Ve = 1-62 X lO'^ cm/s. The BPS result is essentially exact since the plasma coupling is so small. 
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FIG. 16. For the deuterium plasma defined in Fig. 15, the energy (in MeV) as a function of the distance 
traveled (in fiui) is shown for a triton created at threshold with energy Eq = 1.01 MeV in the reaction 
D + D ^ T + p. The solid line is the result from this work (BPS) and the dashed line is the result of Li 
and Petrasso (LP). They give the respective ranges Rbps = 52 fim and i?Lp = 44 /iin, a 20% difference. 
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FIG. 17. Derivative dE{x)/dx of the curve in Fig. 16 split into separate ion (peaked curves) and electron 
components (softly decreasing curves) . This is the triton dE / dx (in MeV / /im) vs. x (in /im) . Again our 
result is the solid curve, and the work of of Li and Petrasso is dashed. The area under each curve gives 
the corresponding energy partition into electrons and ions for this work and that of Li and Petrasso. For 
our result, the energy deposited into electrons is E^^^ = 0.95 MeV and into ions is Ef^^ = 0.063 MeV, 
while LP gives E^^ = 0.93 MeV and ^'j'"^ = 0.075 MeV. These energies sum to the initial triton energy 
Eo = 1.01 MeV. 




F. Plasma Temperature Equilibration 

We now present our results for the temperature equilibration of a plasma in which dif- 
ferent components have different temperatures. This is quite common since plasmas may 
be created in ways that more effectively heat one plasma species over another; for example, 
when a plasma experiences a laser pulse which preferentially heats the light electrons that 
have the larger scattering cross section. We therefore assume that two species a and h are in 
thermal equilibrium with themselves but at two different temperatures and T;,. The rate 
of energy exchange between the subsystems a and h is 

^ = -c.,(r,-r,), (3.57) 

and in this section we present our calculation of the coefficients Cab accurate to order g'^ In Cg^. 

We should remind the reader that there is an hierarchy in which the the light electrons first 
come into equilibrium among themselves, then the heaver ions equilibrate among themselves, 
and lastly the electron and ions equilibrate in temperature. Consider now the typical case of 
a plasma of light electrons and heavy ions, in which the electrons and ions have respectively 
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equilibrated with themselves. The electrons will have a temperature Tg, and and suppose all 
the ions have had time to equilibrate to a common ion temperature Tj. The rate of energy 
exchange between the electrons and the total ion system is 

^ = -Ce:(re-r,) , (3.58) 

where we define Cei = J2i ^ei- With electron and ion specific heats per unit volume Cg and Ci 
defined by dSe = Cg dTg and dSi = CiTi, the rate F is given by 

^ (Te - Ti) = -r (Te - Ti) , (3.59) 

with 

r = C.,(i + i). (3.60) 

In the regime where the electron and ion temperatures are not too different, a case of common 
interest, Eq. (12.12) gives 

TiUle < TeUli : 



where cuf = '^i^f is the sum over all the squared ionic plasma frequencies cuf = e^nj/mj. 

The expression for Cab in the general case is much more complex, and we shall not present 
it here. We will only state that it can be written as a classical contribution plus a quantum 
correction Cab = ^5) + ^^*^' where the quantum piece is given by Eq. (12.49) and the classical 
piece (12.14), C^^ = C^^s + C^b, by the sum of Eqs. (12.20) and (12.26). 



IV. GENERAL FORMULATION 



Thus far we have examined only the de-acceleration of charged particles caused by the 
stopping power dE/dx. This slowing down of a particle keeps it moving in a straight line, and 
all particles starting from the same place with the same velocity slow to a thermal velocity 
at exactly the same final position. But in fact, the width of a narrow beam of particles 
will increase as the particles move through a plasma — there will be a sort of Brownian 
motion in the directions that are transverse to the beam direction. A particle will acquire 
an increasing average squared transverse velocity as it propagates through a plasma. This 
is transverse diffusion. Moreover, the ending positions of a group of particles with identical 
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starting conditions will be spread out along the longitudinal, beam direction to a small 
extent. This is stragghng. 

These random, statistical effects may be accounted for in first approximation by describ- 
ing the charged particle transport by the Fokker- Planck equation 



d 



f{r,P,t), (4.1) 



where /(r, p, t) is the phase space number density for a swarm of particles injected into the 
plasma, each with a common mass m. The momentum derivative d/dp^ acts on everything to 
its right, including the distribution function /. An implicit summation convention in which 
repeated indices arc summed over will be used for vector components, and these indices will 
be denoted by the Latin letters k and /. On the other hand, sums over the plasma species 
h will always be made explicit. In other words, repeated h indices arc not summed, while 
repeated k and / indices arc summed. These conventions are followed in Eq. (4.1). For each 
species index 6, the scattering tensor C^' is symmetric in k and /, as the decomposition (4.18) 
will illustrate. 

The Boltzmann equation reduces to the Fokker-Planck equation when the collisions trans- 
fer only small momenta in comparison with the particle momentum. The way in which this 
limit of the Boltzmann equation works out is described in detail in Appendix^^ C. The 
Lenard-Balescu equation is of the form of a Fokker-Planck equation. Since the right-hand 
side of the Fokker-Planck equation entails an overall derivative in momentum, it conserves 
the particle number density. The terms in the final square brackets on the right-hand side 
of the equation ensure that this side of the equation vanishes for a thermal distribution of 
particles / ~ exp{— /?m v^/2} at inverse temperature provided the background plasma 
components have the common temperature determined by the distribution function, namely 
(3h = {3 for all species h. Thus a thermal distribution of particles is maintained by the 
Fokker-Planck equation. 

Here we shall obtain a precise evaluation of the tensor functions C^' that appear in the 
Fokker-Planck equation. This we shall do to the accuracy of the stopping power dE/dx that 



^In this regard, we should note that although a "diffusion approximation" to the Boltzmann 
equation is outlined in Lifshitz and Pitaevskii [26] , that discussion is in the context of dilute heavy 
particles moving in a gas of light particles. This is a kinematical restriction that is quite different 
than the dynamical case of sharply peaked forward scattering that wc examine in Appendix C. To 
make this distinction clear, we also review the work of Lifshitz and Pitaevskii in Appendix C, where 
we show by a simple example that their result is not internally consistent unless further restrictions 
are imposed. 



37 



we have already discussed. That is, we shall compute not only the leading logarithms, but 
also the constant terms under the logarithms. In this way, among other things, we shall give 
a precise and unambiguous definition of the Landau Collision Integral for Coulomb scatter- 
ing [27]. We shall define the functions C^' by requiring that the Fokker-Planck equation 
reproduce the rate of energy and momentum loss to the plasma species b, quantities that are 
well defined to the leading log plus constant order to which we work, and quantities that we 
compute using the method of dimensional continuation that we have described. 

As a charged particle slows down, large angle scattering events become more important. 
Such hard collisions are not described by the Fokker-Planck equation, and it starts to loose 
its accuracy of describing charged particle trajectories, particularly with regard to the trans- 
verse motion. We shall obtain quantitative criteria for the regions where the Fokker-Planck 
equation ceases to be an accurate description. We shall assess the validity of the Fokker- 
Planck description by computing the rate of transverse energy loss to our order, the order 
which includes the constant terms under the leading Coulomb logarithm. The difference 
between this independently calculated quantity and its evaluation using the Fokker-Planck 
equation tells us when the Fokker-Planck description starts to break down. 



A. Energy and Momentum Transfer Rates 

As we have stated, the coefficients C^' are constrained to produce the energy and momen- 
tum exchange between the charged particle and the background plasma particles of species b. 
To bring this out, we first examine the general case of the transport of some general quantity 
g(p). Averaging over the momentum defines a time dependent spatial density 

QM^ I j^q{p)f{r,p,t), (4.2) 

and fiux vector 

The Fokker-Planck equation (4.1) then expresses {d/dt)Q + 'V ■ !F in terms a momentum 
integral of functions and derivatives acting upon the distribution function /. Hence, by 
partial integration, we may write the result in the form 

|e(M) + V.^(M)^-^J (r,p,<)/(r,p,<), (4.4) 

in which 

C6'(r,P,^)^?(p)- (4.5) 



dQi 
dt 



(r,P,i) = 
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The sign convention has been chosen so that dQb/dt represents the rate at which the quantity 
flows from the projectile to the plasma medium. 

To bring out the meaning of dQb/dt, we note that the Fokker- Planck equation keeps the 
total particle number 



N 

constant in time. Thus 



{,) JV = I d-r j ,(p) /(r, p, (4.7) 

defines a time-dependent, average value of the property q. Wc integrate the local transport 
equation (4.4) over all space. The spatial divergence of the fiux term on the left-hand side 
of Eq. (4.4) is thus removed, and in view of the definitions (4.2) and (4.7) we have 



dt 



-^-E./<'"'/(^^('.P.*)/C-.P.*). (4.8) 



If we now assume that the distribution function / is sharply peaked about a point in phase 
space, at some definite value (rp,Pp), then we may evaluate dQ^/dt at this phase space point 
and take it out of the integral to get 

dt ~ dt ■ ^ ^ 

Therefore, as stated, our sign convention has been chosen so that dQi,/dt is the rate of increase 
of the quantity from the swarm of particles determined by / into the plasma species b. 

Let us now apply these general considerations to the projectile energy and momentum. 
As we shall see, we will get two constraints that completely determine the scattering tensor 
C^''. In this way wc obtain a transport equation that accounts for the secular, long-term 
build up of the changes in the velocity of a charged particle moving in a medium, a transport 
equation with no long or short distance divergences. First we consider the energy density 

f -^-^/(r,p,i), (4.10) 

J (27r)^ 2m ^ ^ ^ ^ 

and energy flux 

S^= f^f'-nr,P,t), (4.11) 



with 



|" + V.5.-E,/^t/(r.P.*). (4.12) 
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According to our general discussion, dE^/dt is the rate of energy loss to plasma species b 
when the charged particle at time t is at the spatial position r with momentum p. Similarly, 
with the momentum density 



and spatial stress 



we have 



T 



kl 



{27ry m 



d 



d^p dPt 



dt' ' ' ' J {2nhY dt 

The general Fokker-Planck evaluation (4.5) therefore gives 

dEb 



f{r,P,t). 



dt 



and 



dPl 
dt 



op'- 



kl 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



Again, the repeated species index b is not summed, while the repeated vector indices k and 
/ are summed over. 

B. Decomposition of the Collision Tensor 

The collision terms in the Boltzmann or Lenard-Balescu equation do not involve gradients 
of spatial variation. Hence the only available vector with which the Fokker-Planck collision 
tensor C^' can be constructed is the particle momentum p = mv. This tensor therefore has 
the general structure^^ 



^^The appearance of the tensor e^'"*{)"* is forbidden by parity invariance. Note that the unit 

vector V, the direction of the velocity, is not wcU-dcfincd in the limit of vanishing velocity, v — 0. 
On the other hand, the tensor C^^ is well defined in this limit of small velocity. Hence a low speed 
constraint must be obeyed: 

v^O : Ab = Bb {l3bv)/2 , 

with 

^;^0 : C^^ = 5^^Bb/2. 
Our evaluations satisfy these constraints. 
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kl 



PbV 2 ^ 



(4.18) 



where the additional factors of l/Pb"^ and 1/2 multiplying the coefficients Ab and Bb have 
been inserted for later convenience. The new scalar coefficients are given by the projections 

1 ^^ .^g^ 



and 



i(i/-l)i3t = C*' 



Alternatively, placing the structure (4.18) in the energy constraint (4.16) produces 



dEb 



d 



V — 



Ab 



(4.20) 



(4.21) 



Likewise, we can relate Bf, to the momentum change dPh/dt in the following manner. Since 
the plasma is isotropic, dPb/dt must point along the direction v. Hence we need only 
compute 



dt 



1 



(4.22) 



where the repeated tensor indices in C" imply a summation over all the spatial axes. In 
view of Eq. (4.16), this can be written as 



.k dPb' 



dEb 



m " dt dt 

The remaining function Bb now can be obtained from Eq. (4.23) together with 



(4.23) 



(4.24) 



The structure of the Fokker- Planck equation guarantees that a swarm of particles with a 
Maxwell-Boltzmann distribution at temperature Tb — 1/Pb remains in thermal equilibrium 
during its interaction with a plasma species b at this same temperature. This aspect can be 
emphasized if we write Eq. (4.21) as 



exp 



We shall find that this formula can be quite convenient for the identification of Ah- Note 
that the thermal average of the rate of energy loss for species b necessarily vanishes. Placing 
the structure (4.25) in the formula for this average gives zero. 
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,2 dEj, 



dt 



0, 



(4.26) 



since the integral entails a total velocity derivative. If there is a common temperature 
T = for all plasma species, then the thermal average of the total rate of energy loss, 



will also vanish. This should be contrasted to the model of Li and Petrasso [24] where the 
thermal average energy exchange (3.30) does not vanish. 



Throughout this work we will often take the distribution function / to be sharply 
peaked in phase space, for example, a distribution peaked about a specific momentum value 
Pp = fripVp. Now, for clarity we write the projectile mass as rup rather than m. Another 
useful case is when / is peaked only about the momentum direction Vp and there is no 
restriction on the absolute value of the momentum itself. For either of these two cases, the 
swarm of particles distributed by / will be called a beam of projectiles, and momentum 
integrals can be performed by the substitutions v ^ Vp and v ^ Vp respectively. Finally, as 
in the preceding example (4.6)-(4.9), we might also consider the distribution function / to 
be sharply peaked in space as well, about a specific point Tp, in which case we can evaluate 
spatial integrals by the substitution r — > Vp. 



Expression (4.21) gives a direct connection between the rate of energy transfer and the 
coefficient Ab, in contrast to Eqs. (4.23) and (4.24) which provide an implicit relation between 
Bi, and the rate of momentum transfer. In the case of a sharply peaked particle beam, 
however, there is also a simple connection between the coefficient and the rate of transverse 
energy flow. As in the second case described in the previous paragraph, we consider a 
distribution / that is sharply peaked in the momentum or velocity direction Vp, and define 
the transverse energy as 




(4.27) 



b 



C. Sharply Peaked Distributions: Projectiles 



1. Transverse Energy 




(4.28) 
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Here our sign convention for this quantity is changed in that q — —E± in Eq. (4.2). Thus, 
the Fokker-Plank evaluation of dE±jj/dt represents the rate of transverse energy flow from 
the plasma to the beam.^''' Equation (4.5) gives 



dE< 



dt 



Since the distribution is sharply peaked about the direction Vp, we can substitute v ^ Vp in 
Eq. (4.29). The first term on the right hand side of Eq. (4.29) involving a single derivative 
vanishes, giving a contribution only from the second term with two derivatives: 



dEj 



dt 



C-{6^-i^^i,l)±=%-2lB,. (4.30) 



F_p 277lp 



Here we have placed an F-P designation on the final result because we will later evaluate 
this transfer rate exactly within our general order of calculation. The comparison of the two 
results provides a signal for the breakdown of the Fokker-Planck equation when larger angle 
collisions become important, coUisions that are not accurately described by the Fokker- 
Planck approximation. As should be expected, we find in Section VI that the difference 
between the Fokker-Planck evaluation (4.30) of the rate of transverse energy transfer and 
the exact rate to our order of accuracy has no large Coulomb logarithm. That is, the 
difference is of relative order of one-over the Coulomb logarithm. We should also hasten to 
mention that in general the transverse spreading of a particle beam is a small effect and so, 
in general, the transverse error is a small error in a small effect. 

2. Velocity Fluctuations 

We have determined the Fokker-Planck coefficient functions C^' by the conditions that 
they correctly describe the rate of energy and momentum transfer between the background 
plasma components b and an arbitrary distribution of test particles or projectiles. An alter- 
native approach was emphasized some time ago by Rosenbluth, MacDonald, and Judd [15]. 
To make contact with this line of development, we consider a distribution function / peaked 



^^To make our sign conventions explicit, we note that we have defined dE/dt and dP/dt as the time 
rate of energy and momentum transferred to the background plasma from the moving projectile. 
Thus, as the projectile slows down, dE/dt and v-dP/dt are positive. On the other hand, Eq. (4.30) 
defines that rate of transverse energy given to the particle by its interactions with the plasma. Thus 
dE±/dt is positive as the projectile slows down. Moreover, Eq. (4.31) describes the change in the 
average velocity of the projectile and so (v) • d{'v)/dt is negative as the projectile slows. 
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at the phase space point (rp, pp), but with finite width. Then, for the averages defined by 
Eq. (4.7), the momentum transfer rate (4.17) imphes that 



di 



rrip dv\, 



(4.31) 



where the scattering tensor C^' is evaluated at (rp,Pp). This resuh requires only that the 
spatial extent of the projectile distribution / is small in comparison with the scale over 
which the plasma properties vary so that C^'' is adequately evaluated at the mean position 
Yp. Since our concern here is with velocity variations, we shall simplify the discussion by 
assuming that the background plasma is spatially uniform so that the spatial coordinate 
dependence of C^' can be entirely neglected. More to the point, the result requires that the 
squared spread in velocity Av'^ is small in comparison with the squared velocity f ^ itself, 
so that C^' may be evaluated at the average (v) = Vp. Here and in subsequent work, we 
neglect the distinction between Vp and (v) except when we are specifically examining the 
velocity fluctuations. [Typically, we will have an equation with {v^ ■ ■ ■) on the left and some 
function F{vp) on the right, by which we implicitly mean F[l(v)).] 
We define 



and next examine 



(4.32) 



(4.33) 



The general relations (4.5) and (4.9) together with the previous result (4.31) and a little 
algebra, show that^® 



dt ' m2 Z-^h ° 



(4.34) 



Writing the Fokker-Planck equation (4.1) as 



d_ 
di 



+ v- V 



demonstrates that it may be expressed as 



(4.35) 



'^^Thc the time rates of change of the average velocity (4.31) and velocity fluctuation (4.34) are 
called "diffusion coefficients" by Spitzer [28]. Our method provides an unambiguous and precise 
evaluation of these to the [In + C] order to which we work. 
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2 dt^ ' 
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Qpk 



m 



4 <^^> 



(4.36) 



This is of the form advocated by Rosenbluth et al. [15], a form also described by Trubnikov 
[16]. These authors work with a completely arbitrary background plasma, a plasma which 
has no aspects of thermal equilibrium. Thus in their work, there is no relationship between 
the rate of change of the average velocity and the squared velocity fluctuation. In our case, 
however, where we assume that the various plasma species h are individually in thermal 
equilibrium at temperature Ti, = 1/Pb, the vector d (f'^) /dt and the tensor d {Av^ Av'') / dt 
are, in view of (4.31) and (4.34), deflned by the same C^' coefficients. When the different 
plasma species have the same temperature T = 1//3, then there is a simple relation between 
the vector and tensor: 



dt^ ' 2 



^ nip dv^p 



^(Av'^Av^) . 
dt ^ ' 



(4.37) 



The coefficients C^' that appear in the Fokker-Planck equation could be determined 
by the contributions of the various plasma species to the rate of velocity fluctuations 
d {Av'^ Av^) / dt rather than by the rate of energy and momentum transfer that we have 
chosen. Such a flxing would give 
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rrir, 
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Pbirip dvl, 



Ab 



(4.38) 



which is precisely our previous determination of the Ab coefficients. The only change with the 
determined by the velocity fluctuations instead of the energy and momentum exchange 
appears in the Bb functions. They would be flxed (in i/ = 3 dimensions) by 



Bb^rup^ 



dEi 
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2 ^^H 
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rrir. 
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V V 
p p 



(4.39) 



As will be shown in Sec. VI, this determination differs from the one that we use by terms that 
are relatively smaller by one over the large Coulomb logarithm. Moreover, the Bb coefficients 
describe only very small corrections to the motion of fast particles. 
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We should stress yet again that we are not working just to the leading order in the large 
Coulomb logarithm, but that we compute exactly the constant terms under this logarithm as 
well: We work to the order g^[ln g'^ + C]. Our method of matching the energy and momentum 
flow determines the coefficients C^' to this order with no ambiguity If instead we would 
have chosen to match to the "diffusion coefficients" d (^Av'^Av^) /dt, then we would obtain, 
to the order to which we work, exactly the same coefficients Ah that determine the rate of 
energy flow, but we would obtain shghtly different coefficients Bi,. The coefficients would 
differ in the constants under the logarithm. These different constants are determined by 
the work of Sec. VI. 

V. TRANSVERSE SPREADING, LONGITUDINAL STRAGGLING 

We turn now to investigate in more detail the nature and effects of the spreading in 
velocity as a projectile moves in the plasma. Although the condition of a small velocity 
spread can be imposed so that it is obeyed initially, it may fail at later times. 

In principle, Eq. (4.31) may be solved to determine (v) as a function of time. With this 
solution inserted in Eq. (4.34) may then, in principle, be integrated to determine the 
fluctuation {v^v'') — {v''){v'') as a function of time. This procedure remains valid so long as 
the fluctuation remains small in comparison with (v)"^. 

The previous decomposition (4.18) evaluated at {v'^) — v^, 

= + (5.1) 

shows that the coefficients Ah determine the longitudinal velocity spread — the 'straggling' 
— while the coefficients B^ describe the transverse velocity spreading, 

and 

Here we should note again that by virtue of the isotropy of the plasma, the average velocity 
(v) = Vp always points along the initial velocity direction Vp (while the magnitude of Vp 
changes with time, its direction remains fixed), and so the average never has a transverse 
component. 

The transverse spreading (5.3) can be expressed as 
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^.V). = (^--^j^y-l^. (5.4) 



dt' ^ P rup dt 

For very fast projectiles, the results (10.43) and (10.45) show that 

m„ „ dE± rrip dE 



^ rUe dt nip dt ^ ^ 



and so 



mpvl » ^ T : {v'^vDt ^ {6'^ - v' v') ^ — {Eo - E) . (5.6) 

^ rUe ^ rup m.p 

Thus at high energies, the transverse angular spreading is of order yjmejvrt^. This is a small 
number for ion projectiles, but of course it is not small for electron projectiles. When an 
ionic projectile slows down to thermal velocities, the transverse velocity fluctuations must 
become of order of the thermal velocity Vp = \/3T/mp. However, until thermal velocities 
are reached, the transverse spreading for ions is always small. 

To assess the nature of the longitudinal fluctuations, the straggling, in a simple way, we 
shall assume that all the plasma species are at a common temperature, (3^^ = T. First we 
examine the motion of a fast projectile so that Eq. (4.38) simplifies to 

dE 

rupvl > T : ^ = Y./"^ ■ (^■'^) 

Thus, since E — mpV^jl^ the rate of longitudinal spreading (5.2) can be written as 

^A»T: I [(„V> 4 1- if , (5.8) 

which integrates to 



mpvl > T : 



[(^ V) - {v') {v')\ ^ = vl *p ^ In (^^) , (5.9) 



where the subscript denotes the initial value. In thermal equilibrium at temperature T, 
the projectile has a root-mean-square velocity Vp — (3T/mp)^/^. Thus the straggling result 
may be written as 

Vp » Vp : [(^V) - {v'^){v)]^ = v'^vl I In , (5.10) 

As the projectile slows to its thermal velocity Vp, its straggling fluctuations become of order 
Vp as are those of all plasma particles. 

The total rate of energy loss is given by 
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The right hand side of Eq. (5.8) can be taken to vanish in the high energy hmit in which terms 
of order T/E can be neglected. Hence with the neglect of terms of this order, Eq's. (5.8) 
and (5.11) imply that^^ (since the direction of (v) is constant in time) 



^ at mr,vi 



dE ^ dEi 



dt dt 



(5.12) 

As we have just discussed, for ions dE±/dt is of relative order uie/mp and can be neglected, 
so that the energy loss rate completely determines the slowing down of the particle. 

Since the logarithm in Eq. (5.10) is a slowly varying function, as the projectile slows 
down from very high velocities to speeds that are more nearly of the order of the thermal 
velocity, it acquires velocity fluctuations in the longitudinal direction that arc only slightly 
larger than Vp. This justifies the integration of the slowing down equations (5.11) and (5.12) 
from very high velocities to just above thermal speed using the average velocity (f ) = Vp, 
since the velocity spreading is relatively very small. However, this simple picture of the 
essentially deterministic motion of an individual particle breaks down when the particle 
speed approaches the thermal speed of a particle of its mass. In this region, a statistical 
distribution of particles must be employed as the proper description. The Fokker-Planck 
equation can be used to describe the time evolution of the phase space density /(r, p, t), and 
our computation of the coefficients C^' that enter into the Fokker-Planck equation remain 
valid. What breaks down is the single-particle description of a particle losing well defined 
amounts of energy and momentum. 

What we have just said means that whenever l/PhnipVp corrections become important, 
the notion of a well-defined projectile trajectory breaks down. Hence the corrections given 
in the second part of Eq. (4.21) are never relevant for the description of a single particle 
which, in any relevant region, is described by 

dE 

rUpvl > T : ^ = VpAb , (5.13) 

which also determines the complete motion of the particle. Nonetheless, we have used the 
form (5.13) in all of our calculations of dE/dx = {1/vp) dE/dt. 

We have noted that the transverse spreading of an electron projectile may be significant, 
even for fast particles. To emphasize this point, we quote the high speed limits for the total 
and perpendicular energy loss of an electron that follow from Eq's. (10.43) and (10.45): 



^^This result also follows from Eq. (4.37) as it must. 
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mevl:$>T: ^^in , (5.14) 



dt 47r I'p \ tiWe 



and 



rUeV^ > T : — r— = > — - — - — In — . (5.15) 

Here is the electron plasma frequency defined by = e^ne/rUe while — elrih/Th. The 
general case is sufficiently well illustrated by the specific case of a fully ionized hydrogen 
plasma with equal numbers of electrons and protons, rig = Up, and equal proton-electron 
temperatures, = Tp = T. The reduced masses are given by rriee — 'me/ 2 and, with the 
neglect of the small electron proton mass ratio, rUep — mg. Hence 



meV„ > T 



dE^ ujj , / 2mlvl 
= ^ In ' 



dt 47r Vp \h?l3e'^2n 



ujI [mevl T 



In ^ . (5.16) 

47r Vp \ nujf. nwej 

For the plasma parameters that we have used above, \n{T/tujJe) is not a large (or small) 
number. Hence dE±/dt is about the same size as dE/dt, and so the the electrons do not 
slow down along a straight line. The electron motion in a plasma requires the use of a 
Fokker- Planck description of an ensemble of particles. 

VI. VALIDITY RANGE OF THE TRANSPORT EQUATION 

The exact — to our order — rate at which the transverse energy of a projectile increases 
is the sum of the leading z/ < 3 result computed from the Lenard-Balescu equation plus the 
leading z/ > 3 result computed from the Boltzmann equation. The difference of this with the 
evaluation (4.30) given by the Fokker-Planck equation provides a signal for the breakdown 
of the Fokker-Planck description. This measure is 



dE\ 



dt 



dEi 



exact 



(6.1) 

F-P 



It is worthwhile providing here the results of this assessment, detailed in Sec. XI below, 
so as to conclude our general review in a unified manner. But before presenting these results, 
some general remarks may help clarify what we are doing. We have used the two functions 
Y-dPb/dt and dEb/dt as inputs to determine the two scalar coefficients At and Bb that define 
the transport tensor C^' that appears in the Fokker-Planck equation. Then the time rate of 



49 



change of the perpendicular energy, dE±/dt, may be found from C^', a determination that 
we denote by the F — P label. The point is that, within the Fokker-Planck approximation, 
only two of the three functions v • dPb/dt, dEb/dt, and dEx/dt are independent functions. 
However, if the Fokker-Planck approximation is not made, then these three functions are 
linearly independent functions. The difference (6.1) is thus a measure of the error in the 
Fokker-Planck description. 

Since the p < 3 contribution to the Fokker-Planck coefficient C^' is the same Lenard- 
Balescu equation that is used to evaluate the transverse energy in this region, the difference 
defining is given by just the 1/ > 3 parts, 

dE> 



dt 



exact 



dt 



(6.2) 



F-P 



with both terms computed from the scattering cross section formula that is equivalent to 
the Boltzmann equation as is described in Sec. VIII. This computation is given in detail in 
Sec. XI, with the result (11.21) that 



4:71 2mr. 



1 



(6.3) 



In the high speed limit, ■mbVp/2 ^ T^, the exponential is highly damped. Hence, in this 
limit, the upper integration limit may be extended to — > 00, and one finds that 



00 



(6.4) 



47r 2(3bmpVp 

No "Coulomb logarithm" appears here because the difference A^ is not sensitive to small 
angle scattering or, equivalently, to large distance collisions. 

In the low speed limit, the exponential may be expanded in powers of v'^. The integrations 
involved in the zeroth order term vanish, and one finds that 
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47r 15m, 



1/2 



PbUlbV, 



The definition of E±b involves the tensor VpVp, which is undefined as Vp 
must be accompanied an additional factor of Vp (as VpVp = Vp), thereby giving a we 
and quadratically vanishing tensor VpVp. 

These hmits can be used to assess the validity of the Fokker-Planck evaluation 



(6.5) 



0. Therefore it 
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(6.6) 
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The high- velocity hmit (10.45) for St gives 



oo 



dEi 



dt 



F-P 



47r PbTUpVp \ Hko 



(6.7) 



As should have been expected, this is larger than the difference A;, by the logarithmic factor 



2 1 /2m^\ _ 



As far as the transverse spreading is concerned, this shows, at least at high energy, that 
the Fokker-Planck description is valid only to leading logarithmic order: The validity of the 
transverse spreading given by the Fokker-Planck equation is valid to the accuracy to which 
the Coulomb logarithm is large in comparison to unity. It must immediately be remarked, 
however, that, in the high-energy limit, the rate of energy loss 



oo 



dE_ 
It 



47r PerrieVj: 



In 



(6.8) 



is much larger that the rate of transverse spreading, larger by the very large ion/electron 
mass ratio nib/me. Thus the spreading entails very small angles, and one can tolerate a 
rather large error in this small effect. 

In general, the error is smaller than the transverse energy spread itself by a factor 
of one over the Coulomb logarithm. To see this, we may, for example, turn to Eq. (9.5). 
It contains a logarithm whose factor K in its argument cancels against that in another 
contribution to yield the large Coulomb logarithm [with a quantum or classical cutoff as is 
appropriate to the velocity Vp] Thus the leading Coulomb logarithm contribution to the 
Fokker-Planck transverse energy rate (6.6) is given by 

1/2 



dt 



An \2nPb 



rrir 



J duu (1 — m) exp | — 



'brribVp u 



(6.9) 



F-Plog 

in which Lb is the appropriate Coulomb logarithm for this plasma species h. This is indeed 
larger than A5 [Eq. (6.3)] by essentially the factor Lb- 

The result (9.12) gives the small velocity limit of the Fokker-Planck approximation (6.6): 
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(6.10) 
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This result is in accord with the comments of footnote 16, which explains why 



Vp^Q : CI' const • 5'^' = /2 . (6.11) 

Using Eqs. (4.16) and (6.11), along with Eq. (4.30), in an arbitrary number of spatial di- 
mensions V we find 



Vp^^ 



dt 2mp ' dt 2mp ' 



and hence 



dE^b v-l dEb 

Vj,^ u : — - — = ;— . (6.13) 

^ dt V dt ^ ' 

Our sign conventions, in which dEb/dt is the rate of energy loss of plasma species 6, while 
dEj^b is rate of transverse energy gain of species 6, dictate the relative minus sign above. 
The low speed limit of dEb/dt for = 3, as computed by Eq's. (4.21) and (9.9), is indeed 
just the factor (3/2) times the result (6.10) for dE±b/dt. On other other hand, the result 
(6.10) is quite different than the low speed hmit (6.5) of the error which behaves as Vp 
when Vp 0, not as a constant as given in Eq. (6.10). Thus, the Fokker-Planck equation 
gives a very accurate description of the transverse spreading of low velocity particles. 
We turn now to the details of our calculation. 

VII. LONG DISTANCE EFFECTS DOMINATE WHEN u <3 

When the spatial dimensions are less than three, long-distance, collective effects are 
dominant. This "soft physics" is described to leading order in the plasma density by the 
Lenard [6]-Balescu [7] equation. ^° Indeed, to leading order in the density, one can prove that 
the rigorous BBGKY hierarchy reduces to the Lenard-Balescu equation when u <3. This 
can be demonstrated, for example, by carefully examining the discussion given in Nicholson 
[20] or in Clemmow and Dougherty [10]. It is significant that the proof of the reduction 
of the BBGKY hierarchy to the Lenard-Balescu equation breaks down at precisely u — 3: 
this happens because of the appearance of short-distance, ultra-violet divergences, which are 
absent in dimensions less than three. 



^° Again we note that Refs. [10,19,20] contain well written expositions. 
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The Lenard-Balescu equation for the case of interest in which each background plasma 
species b is in thermal equilibrium and described by a Maxwell-Boltzmann distribution at 
temperature T], = l/Z^b is of the Fokker-Planck^^ form (4.1), using 

with^^ Vp = Pp/rrip and v;, = pi^/nii, the velocities of the projectile and of the background 
plasma species b. The collective behavior of the plasma enters through its dielectric function 
e(k, u;). For a dilute plasma, the case to which the Lenard-Balescu equation applies, the 
dielectric function is given by^^ 

where the prescription 77 — > 0+ is implicit and defines the correct retarded response. 

A. Projectile Motion in an Equilibrium Plasma 

Taking the projection (4.19) of Eq. (7.1) and setting k • Vp = kcos9, or alternatively 
taking the trace of the tensor indices in Eq. (7.1), yields 



2 2 



(2^^|e(k,vp-k) 



(7.3) 



Here the less-than superscripts on and C"^ are written to make it explicit that we 
are now working in spatial dimensions strictly less than three. Using the delta function in 
Eq. (7.3) to remove the component of Pb = m^Vf, along the k direction, and then integrating 
out the remaining 1/ — 1 components of using the Maxwell-Boltzmann distribution 



^■•^ Although this is a purely classical result, the factor of h~'^ in the measure (F-pb/ {2'KhY is used to 
convert the momentum integral of the dimensionless phase-space distribution fb{Pb) into a particle 
number density. 

^^We now use the subscript p to distinguish the projectile velocity Vp (which we previously simply 
denoted by the unadorned v) from the velocities of the background plasma particles of species 
b. 

^^See, for example. Section 29 of Ref. [21]. 
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JfelPfe) = exp < -— rribVf, } (7.4) 



\ rrib J 12 



reduces Eq. (7.3) to 



= 4 



exp < - — mfe f „ cos^ 9 > ^ { cos^ ^ , 1 ^ , 



where 



(27r)- V 27r/56 " 1^ 2 " ^ J |A;2e(k, A;^;^ cos ( 

(7.5) 



= /?6 rift (7.6) 



is the contribution of species h to the squared Debye wave number. 

To work out this result, we first note that the structure of the dielectric function (7.2) 
can be simplified. We use the explicit Maxwell-Boltzmann form for the distribution func- 
tion /c(Pc) to compute the derivative in Eq. (7.2) and then integrate out the momentum 
components of Pc that are perpendicular to k. This gives the structure 

k'^ e(k, kvp cos 9) ^ k^ + F(vp cos 9) . (7.7) 

The F function appears in the form of a dispersion relation 

with the spectral weight 

Ptotal(^^) = Y] Pc{v) , (7.9) 
'c 

where 



Pc{v) = Klv^(^^exp}^~l3cmcv'^j. (7.10) 

For future use, we note that F satisfies the relations 

F{-u) = F*{u) (7.11) 

and 

Im F{u) = i [F{u) - F*{u)] = 7rptotai(^) ■ (7.12) 
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As a first application of this structure of the dielectric function, we note that for a plasma 
all of whose components have the same temperature, the total energy loss reads 



where 



dE< 
dx 
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dx 



\_ T d 
rripV dv^ 



(7.13) 



(7.14) 
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we have 



27r 
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|/c^e(k, kvp cos 9) 
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Im F{vp cos 9) 
\k'^ + F(vp cos 9)1^ 



— Ini 

k e{k,Vpk cos 9) ' 



- — — — cos 9 Im — ; — 

{27:)" k e{k,Vpk cos 9) 



(7.15) 



(7.16) 



Except for the term involving the derivative in the energy loss formula (7.13), Eq. (7.16) is 
just the energy loss to Joule heating the plasma, the energy loss obtained by using Fourier 
transform techniques to compute the volume integral of j • E, where j is the current of a point 
particle moving with velocity Vp and E is the electric field produced by this current. The 
additional term involving the derivative provided by the correct Lenard-Balescu transport 
equation ensures that the total energy loss vanishes for a swarm of particles with a thermal 
distribution of velocities at temperature T — 1//3. 



B. Calculating the Coefficients 

We use Eq. (7.7) and the results that follow it to place the energy loss coefficients (7.5) 
in the form 

(.<J_ nii<\=il f kpb{vpCos9) 
X'PbVp' ' j PbVp J {2ny \k'^ + F(vpC0s9)f 

The wave number integration may be performed by passing to hypcr-spherical coordinates. 
For functions depending only upon the radial coordinate k and the polar angle 9, we may 
write 



jcos^,^), (7.17) 
I cos9 \ 
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where il,^_2 is the sohd angle subtended by a, {u — 2)-dimensional sphere. The /c-integral 
in Eq. (7.5) is of the form 

fOO 7,1/ 

which is finite for z/ < 3 and log-divergent at z/ = 3. Despite the fact that one thinks in terms 
of integer dimensions, one is nonetheless free to perform the integral (7.19) treating u as an 
arbitrary complex number. Moreover, the solid angle factor fij,_2 has an analytic form that 
extends to arbitrary complex dimensionality u, and the power sin'^"^ in the polar angular 
integration can also obviously be extended to arbitrary complex u. The whole expression for 
the energy loss rate can be extended to a space of arbitrary complex dimensionality u. The 
physical dimension z/ = 3 is, however, a singular point, namely a simple pole. Nonetheless, 
we can regularize this infinity (that is, render it finite) by formally treating u as complex 
number differing slightly from three. In any well-defined physical process, all terms that 
diverge in the i/ 3 limit must cancel among themselves. For the problem at hand, we will 
show in the next section that short-distance scattering, which has not yet been included, 
produces a divergence as z/ — > 3 that exactly cancels the aforementioned divergence. This 
renders the experimentally measurable energy and momentum loss finite in three dimensions. 

Let us now evaluate the integral (7.19). It is convenient to add and subtract a (well 
chosen) term so as to express the integral as a sum of two pieces, the first having a divergence 
when z/ ^ 3 but with no 6 dependence, the second a finite term which does have the rather 
complicated ^-dependence of the function F{vpC0s9): 

1 1 

{k^ + F){k^ + F*) ~ k^ + K^\ ' 

(7.20) 

where K is an arbitrary (^-independent) wave number. Since the final result cannot depend 
upon K (as we have merely added and subtracted the same A'-dependent quantity), we can 
choose its value as a matter of convenience. The first term of Eq. (7.20) is 6'-independent 



^0 



dk 



\k^ + F\ 
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dk 



k" 
k^ + K^ 



+ 



r dk 



In general, by a "d-dimensional sphere" we mean a sphere whose hyper-surface is of dimension 
d, which can be thought of as a sphere embedded in (d + l)-dimensional Euclidean space. Points 
on such a sphere centered at the origin with unit radius satisfy Yle=i ~ ^- '^^^ solid angle 0^ is 
simply the surface area of this unit sphere, and it can be expressed as = 27r('^+i)/Vr((d+ 1)/2). 
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and divergent as — > 3; the second term is finite in this hmit (so we have taken u — 3), but 
its 9 dependence is non-trivial. 

The first integral on right hand side of Eq. (7.20) is straightforward to evaluate^^ 

1 K"-^ 



dk 



K 



1/-3 



TT 



+ 0(1/ -3). 



(7.21) 



k^ + K^ 4 sin (tt^) 3 - i/ 

By partial fractions, the second integral on the right-hand side of Eq. (7.20) is easily evaluated 
in terms of logarithms. Thus, in the z/ 3 limit. 
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(7.22) 



|fc2 + F|2 3-1/ 2{F-F* 

The derivative of Eq. (7.22) with respect to K vanishes when z/ 3, hence the X-dependence 
cancels in the physical limit. With the aid of Eq. (7.22), the coefficient functions now appear 
as 
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PbVp 
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PbVp (27r) 



^^-^ / (icos6' sin'^~'^6'/0fe(f„cos/ 
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Fin ( 



cos 9 , 



cos^ 



(7.23) 



where Eq's. (7.10) and (7.12) have been used to simplify the notation. 

The second set of terms in the curly braces in Eq. (7.23) are obviously finite in the i/ — > 3 
limit. We take this limit for these terms, and write a decomposition into singular and regular 
parts. 



The singular part is given by 

_ el fKy-^ 1 1 



— *^6,s + 



^6,R ■ 



(7.24) 
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47r 2n \27r J 3 — u jS^Vp J_i 



+1 



d cos 9 sin^ ^ 9 Pb{vp cos 9) ^ cos 9 , 



cos^ 



(7.25) 



The regular part can be simplified slightly by using the refiection property F[—u) — F[uy 
noted in Eq. (7.11) while the ratio Paiu) / ptota\{u) is even in u. Hence 



2^ One sets /c^ = x K'^ and expresses the resulting integral in terms of a contour integral involving 
the discontinuity of x^'^~^^^^. The contour integral may then be opened up to enclose only the 
simple pole at a:; = — 1, which gives the result (7.21). 
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47r /SftVp 27r 7_i Ptotai('ypcos6') V -^^ 7 1 ' cos6' 

(7.26) 

The corresponding coefficient functions summed over all the species in the plasma are easily 
obtained for a plasma at a common temperature, since one only needs to place 

T = 1 (7.27) 

in Eq. (7.26). To write the singular contribution (7.25) in a convenient form to combine with 
the u > 3 result, we change variables to u — cos^ 9 to get 



C. Asymptotic Results for Large and Small Velocity 

The asymptotic forms of our results for large and small projectile velocities Vp are of 
interest. Here we shall work out these limits for the regular terms. The corresponding limits 
for the singular terms are much easier to compute once they arc combined with the singular 
terms produced from the z/ > 3 calculation, and so we defer this until later on. 

To obtain the small velocity behavior of the dielectric function, we first add and subtract 
VpCosO/v in the numerator of the integrand of (7.8) to get 

F(^;.cos^) = 4-y kI r^dv "^^^ ./^exp|--/5emyl , (7.29) 

' " " Vpcos0 + ir]-v\l 27r ^\ 2^' ' J ' ^ ' 



where 



4 = E '^c (7.30) 



is the total squared Debye wave number of the plasma. We now make use of the relation 

^ — ^ -iT^S{vpCose -v)+V \ , (7.31) 



Vp COS 6 — V -\- IT] ^ Vp COS 9 — V ^ 
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in which V denotes the principal part prescription. Since V{l/x) defines an odd function, 
the translation u — v — Vp cos 9 of the integration variable gives^^ 



F{vp cos I 



cos I 



exp I -^Pcruc {vp cos 61) ^ I 



/ 

Jo 



sinh {Pc'iTT'c u '^p cos 9) exp < —-(3cmc u 



(7.32) 



10 u 

In this form the small Vp limit is reduced to the evaluation of elementary Gaussian integrals 
and we have 



: 



F{vp cos 9) ^ K^-^ kI (5c rric cos^ 9 + 0{Vp) + m ptot3x{vp cos 9) , 



(7.33) 



where we note that ptotaK^*) starts out at order u. Placing this result in Eq. (7.26) produces 



and 



Vp^O 



b.R 



e'p 2 fPb^bY' f 



In 



TT v:, 
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^1+1 
2 



(7.34) 



1/2 



1 - -^Pbmbvl 



In 



4/5c-c + ^ [E 4(/^cm.)H'|. (7.35) 



26i 



The integral defining the real part which appears here may be written in the form 

+ 00 



„ , 1 /■+~ dv , 
'^y^j ~ o — smnzra 



since the integrand in an even function. Differentiating this with respect to x, writing out the 
resulting hyperboHc cosine in exponential terms, and completing the square yields two simple 
Gaussian integrals which give d'{x) = y^expjx^}. Thus we have the alternative evaluation 

d{x) = Vtt / dy exp{y^} , 

^0 

which is essentially Dawson's integral. 
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To obtain the large projectile velocity limit of the ^^j^ coefficient, we first note 
that the numerator and each term in the denominator of the spectral weight ratio 
pi,{vp cos 9) / ptotai{vp cos 9) [where pj, and ptotai are defined in Eq's. (7.10) and (7.9)] con- 
tains a factor exp{—(3cmcVp cos^ 9}. In view of the very small electron/ion mass ratio, these 
exponential factors approach much faster for ions than for electrons in the Vp ^ oo limit. 
Thus the spectral weight ratio is very small except for the case in which the index b refers 
to the electron, b — e, for which case we have 

PeivpC0s9) 

Vp^oo : f — 1 . 7.36) 

Ptota\{VpCOs9) 

Thus the sum defining the energy loss to all plasma particle species is dominated by the 
electron contribution in the large projectile velocity limit: 



Vp ^ oo 



At = yAt^^A<. (7.37) 



The asymptotic limit that we are about to obtain is valid when f3emeVp/2 1, or when 



1 Tfl 

Ep= -mpv^:^ ^T,. (7.38) 

^ TTtp 



To obtain this large velocity limit, we first replace the spectral weight ratio by unity in 
formula (7.26). We then write cos^ = 2; in Eq. (7.26) and note that since the integrand 
is analytic in the upper-half z plane, we may deform the original integration along the 
— 1 < 2; < -|-1 portion of the real axis into a semicircle of unit radius in the upper-half z 
plane. Thus (7.26) becomes 

When M = is a large complex variable in Eq. (7.8) we can perform simple Gaussian integrals 
to obtain the limit 



This limit actually entails a sum of terms over all the species 6, with the electron mass irie 
replaced by that of the species 6, mg mi,. However, because of the very large ratio of the 
ion masses to the electron mass, this sum is dominated by the electron contribution that we 
have written. Keeping only the dominant electron terms and placing this limiting behavior 
in Eq. (7.39) gives 
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i'p — > oo 



47r 2(3emevl \ nl 

The asymptotic behavior of C"^ is quite different. The final factor of 1/ cos 9 in Eq. (7.26) 
emphasizes the region about cos^ = 0. Hence the large Vp considerations given above for 
A^^ do not hold. Instead, we must regulate the integrand near cos^ = 0. Since F(0) = k^, 
this is done by writing Eq. (7.26) as 

(jiiK^^J: L r\coce ^^^^''^^^^^ F{vp COS 9) 

'''^ AtT PbVp 2% J_i cos 9 ptotal (Vp COS 9) 

For the integral involving the first constant logarithm, we note that the overall function 
multiplying F(t'pCos6') is odd in cos 6'. Hence in view of Eq. (7.12), we may replace 

F{vp cos [Fivp cos 9) - F{-Vp cos 9)] = 7riptotai(^^p cos 9) (7.43) 

in this part. Since In {F{vpC0s9)/ F{0)) vanishes at cos^^ = 0, the integral involving it may 
be treated in the same way as was done previously for A^,^. Only the electrons in the plasma 
contribute to this second piece. Thus, using the Kroenecker delta function S^^e to distinguish 
the electron contribution, the two parts give 

ii<_ 4 ^ ^ r\^^.nPbi'"pC0s9) f kI 



C R = -iri^o / c^cos^ ""^ . 1^ \ 172 
A7r(3bVp2j_i cos9 

^^.4^rS^(«,e.*)l.(£^). (7.44) 



AnpbVpJo 27r ^ ^ ^ V ^(0) 

In view of the exponential damping in the definition (7.10), we may replace the cos^ inte- 
gration limits of ±1 by ±oo to obtain the asymptotic form 



Vp ^ oo : 



j dcos9 — > J dzK,t,Vp^ exp <~-Pbmi,VpZ | = . (7.45) 



On the other hand, in the large velocity limit. 



Vp ^ 00 : 



r^Five^'^) In(^^fe^) ^ (-^e--A In 
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This is of relative order Te/nieVp and hence may be neglected in the asymptotic limit, leaving 
only 



oo 



r<ii< _ 




(7.47) 



VIII. SHORT DISTANCE EFFECTS DOMINANT WHEN ly > 3: 

CLASSICAL CASE 



To the leading order in the plasma density with which we are concerned, the Boltzmann 
equation correctly describes the Coulomb interactions in the plasma for spatial dimensions 
u larger than three. Again one can prove that, to leading order in the plasma density, the 
rigorous BBGKY hierarchy reduces to the Boltzmann equation when the spatial dimension 
exceeds three. The Boltzmann equation for the phase-space density fa{Pa) of species a 
reads 



We suppress the common space and time coordinates r, t and write the collision term involv- 
ing species b in the form 



Here, although at this stage the scattering process is taken to be purely classical, a quantum- 
mechanical notation has been adopted^^ to describe the scattering of the particles of mass 
rUa and rrib, the scattering from the initial momenta Pa = rUaVa, Pb = "m-b^b to the final 



^'^This can be established by carefully examining, for example, the discussion of the derivation of 
the Boltzmann equation from the BBGKY hierarchy given in Sec. 3.5 of Huang [18]. The derivation 
breaks down at = 3 because of the long range of the Coulomb force. 

^^The roles of the factors of Planck's constant h that appear in Eq. (8.2) are worth pointing 
out. It suffices to consider the factors associated with the first product of phase-space densities in 
the square brackets in Eq. (8.2). The two factors of h~'^ that appear in the first two integration 




(8.1) 



Cab{Pa) = 




fb{p'b)fa{p'a) - MPb)fa{Pa) ■ 



(8.2) 
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momenta p'„ = fnay'a, p't = "^feV^, with the scattering amphtude T{W, q^) depending on 
the center-of-mass energy W and the squared momentum transfer q^. It is convenient to 
employ this quantum-mechanical notation even for classical scattering for several reasons. 
It explicitly displays the complete kinematical character of a scattering process, including 
the detailed balance symmetry. It explicitly shows that the Boltzmann equation may be 
generalized to an arbitrary number of spatial dimensions v. It connects the collision term 
explicitly with the cross section generalized to v dimensions as shown in Eq. (8.8) below. 
Finally, it shows that the collision term (8.2) vanishes when all the particles are in thermal 
equilibrium with the generic densities /(p) ~ exp{— /?mv^/2} because of the conservation of 
energy enforced by the delta function. 

A. Projectile Motion in an Equilibrium Plasma 

The transport of energy, momentum, and transverse energy was discussed at the start 
of Sec. IV. To do this in terms of the Boltzmann equation, wc first review the standard 
treatment for the sake of clarity and completeness, and to establish our notation. In general, 
we deal with a momentum-dependent quantity g(p) which gives a spatial density for species 
a as 

Qa(r, t)^ j g(p„) /„(r, p„, t) , (8.3) 

and flux vector 



volume elements are the factors of h that appear even in purely classical statistical mechanics. They 
change the dimension of d'^p' to that of an inverse volume (in a i/-dimensional space) so that the 
momentum integral of /(p') becomes a particle number density. The remaining factor of h"" in 
the d'^Pb measure just cancels the conventional factor of W associated with the total-momentum- 
conserving delta function. Since the dimension of d{x) is x~^, the single factor of h associated with 
the energy-conserving delta function produces a quantity with the dimensions of time. So far, we 
have the dimension count L^^^ T. The final factor to be examined is the scattering amplitude. To 
obtain its dimensions and its overall h dependence, we consider the first Born approximation result 
T = V{q/h), where V is the Fourier transform of the potential and q is the momentum transfer 
in the scattering. Thus the scattering amplitude has the dimensions L'^, and we conclude that 
the collision term Cab{Pa) has the dimensions of a rate, T"^, as it must. As we shall later see 
explicitly, all the factors of Planck's constant h cancel in the classical limit save for those that 
convert the two momentum integrals of the phase-space densities into number densities, the h 
factors that appear even in classical statistical mechanics. 
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The Boltzmann equation (8.1) then gives 

^^Qa{r,t) + V':F'^{r,t) = I ^^QiPa)J2,^a,{r,Pa,t). (8.5) 

The sum involves colhsions with the other particle species in the plasma. Thus, we have 
an unambiguous identification of the rate of transfer to each species b in the plasma. The 
scattering amplitude, delta functions, and momentum integrations in the collision term are 
symmetrical under the interchange of initial and final particles. Hence, we may make the 
replacement 

?(Pa) fbiPb)faip'a) - MPb)faiPa) ^ qip'a) - QiPa)] fb{Pb) faiPa) ■ (8.6) 

The collision term in Eq. (8.5) can be expressed as 

/ (^ ^b ^^'^'^ ^) = / ^ > (8-7) 

where the sign is chosen so that dQab/dt gives the rate at which the property Q is transfered 
from species a to the plasma species b. We now concentrate on a particular momentum 
and position of a particular "projectile" particle p and identify dQb/dt as the transfer from 
this projectile particle to the plasma species b. Now, in general, the cross section for the 
scattering of particles p and b into a restricted momentum interval A is given by 

(27rn) S (^rUpVp^ + ^^m^v'^ - ^rUpVp^ - ^mbVb^^ . (8.8) 

Using this definition, we find that 
dQt _ f d'^pb 



dt 

Here 



/d'^Pb f 
J^:^fbiPb) Vpb J dapb [q{Pp) - qipp)] . (8.9) 



Vp6 = Vp - Vfe (8.10) 



is the relative velocity between the incident particle and plasma species b, with magnitude 
"^pfe = Wpb\- Henceforth we use the > superscript to emphasize that this is the leading result 
for z/ > 3. Except that we work in a space of arbitrary dimensionality u, the result (8.9) is 
the familiar one: The rate of change of a quantity is its change in a collision times the rate 
of these collisions — which is given by the cross section folded over the incident fiux values. 
In this form, rather than the change of the quantity brought about by the scattering into 
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and out of a momentum region as described by the Boltzmann collision term, a relabeling 
of variables expresses the rate of change in terms of the change in each coUision. The result 
(8.9) expresses the rate in an obvious form, but it does not make manifest the fact that the 
rate (8.9) for a quantity that is conserved in the collision vanishes when integrated over a 
thermal, Boltzmann distribution with the same temperature of the plasma. 

Some momentum integrations may be performed by passing to the center-of-mass coor- 
dinates, where the total and relative momenta are defined by 

P = mpVp + mfeV6, (8.11) 

p = nipb-Vpb = -T^ {mbPp - rripPb) , (8.12) 
^vipb 

with Mpb the total mass and rripb the reduced masses of the system, 

111 

Mpb ^m.p + mb, = 1 . (8.13) 

nipb nip nib 

Similar expressions hold for the final state variables. We can now write the momentum and 
energy conserving delta-functions as 



(1111 
-mpv'J" + -nibv'^^ - -nipVp^ - - 



nibVb 



= b^""^ (P' - P) 5 f — ^1 , (8.14) 

^ ' \2nipb 2mpb) ^ ' 

and since there is a unit Jacobian in passing to center-of-mass coordinates, d^Pb d^v'p = 
d'^P'd'^p', we have 



/ dapb= I \T{W,q^)? {27rh)6(^ ^] . (8.15) 

7a yA(27r^)"' ^ '^^1 ^ ' \2nipb 2nipb) ^ ' 



V. 



We now note the energy in the center of mass is given by 

W^^, (8.16) 
2mpb 

while the momentum transfer in the scattering is 

q = Pp - Pp = Pb - Pfe = p' - P ■ (8.17) 

Let us first apply these considerations to the rate of energy transfer dEb/dt. This is 
obtained from the general formula (8.9) with 

q{Pp) - q{Pp) ^ Irnp [V^' - vj] . (8.18) 
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The conservation of momentum P' = P, and the energy constraint p'^ = p^, allows us to 
write the energy change of the projectile as 

im,«-.^)=^P.q. (8.19) 

Since the scattering in the center of mass frame is axially symmetric about the initial mo- 
mentum p, the transverse components of q average to zero in the scattering process, and so 
we may make the replacement 

P-q^4(P-P) (P-q)^-AP-P?', (8.20) 

with the last form following from the energy constraint p'^ = p^. Thus we may write the 
energy loss as 

dE> _ 1 f d'^p, „ , , P-p f, 2 



dx Vp 



To extract from Eq. (8.21), we use the relation (4.25), which we repeat here for 
convenience: 

exp Pbrripvl^ PbrripVp ^ ^ *p Pbrripvl^ A • (8.22) 

Since P = m^Vp + nibYb and p = rUpbVpb, and since the cross section integral is only a 
function of the magnitude of the relative velocity Vpb, the integrand of Eq. (8.21), multiplied 
by the exponential factor in this relation, has the velocity dependence of the form 

-Pb exp |-i Pbmpvl^ fb{pb) P • pX{vpb) 

~ X{vpb) rUpb-Vpb ■ + cxp (5b {mpvl + rribvl) | . (8.23) 

Integration by parts in the pb integral replaces the action of d/dvb on the exponential factor 
by —{d/d'Vb) • WpbX{vpb). Since = Vp — Vf,, this action of —{d/dvb) is equivalent to that 
of {d/dvp), and so we have, effectively within the p^ integral, 

exp I - ^ Pbrripvl | (pb) P • p X {vpb) • exp | - ^ (^brripvl | /(, (p^) p X {vpb) . 

(8.24) 

The only direction produced by the p;, integral is that along Vp. Hence we may replace 
P ^ Vp (vp • p). We thus arrive at the structure (8.22) upon the identification 
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1 f d'^Pb f 



2mb J (271/1) 
Here we have simplified an overall factor by using 



- ' (8.26) 



Mpb rupb rub ' 

The remaining independent coefiicient may be taken to be C^' which, according to 
Eq. (4.23), is given for > 3, by 

^ r<ll> — .,fc "-^6 "-^b ('8 97^ 



This difference involves 

Vp • mp (v; - Vp) - ^ mp (^v;' - vfj = -^m,, (v; - v^)' = -^^^ ■ (8-28) 
Thus we rather quickly find that 

= ^ / A(Pb) / da,b . (8.29) 
In summary, we may write the results as 



B. Classical Coulomb Scattering 

For the remainder of this Section we shall treat only the case of classical scattering. We 
defer the discussion of quantum-mechanical corrections to Section 10. 

1. Cross Section Integral 

We now apply the energy loss formula to the case of classical Coulomb scattering which 
results from the reduction of the classical BBGKY hierarchy. In v dimensions, the element 
of differential classical cross section is given by 

da'p^ = n,_2h''-^dh, (8.31) 

where is the area of the unit u — 2 sphere and h is the classical impact parameter. 

Hence 
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oo 

da^, = fi,_2 j db b'^-y (b) , (8.32) 

with the momentum transfer related to the scattering angle by 

q\b) =4p^sm^e{b)/2) , (8.33) 

where 9{b) is the scattering angle as a function of 6 in dimensions. The classical planar tra- 
jectory of a particle moving in a central potential is independent of the spatial dimensionality 
u. Thus the familiar formula for the scattering angle 

-1/2 

(8.34) 



holds for arbitrary spatial dimensionality u. Here rmin is the lower turning radius at which 
the angular brackets in the integrand vanish. The Coulomb potential energy in u dimensions 
is given by 

V(r)^.,^.(r) = ;^;;^;^>, ^. (8.35) 

which follows from Gauss's law in i/ dimensions, 

e, = [d^rV-B,^ -Q,_i r"-' . (8.36) 

J dr 

Changing the integration variable to m = 1/r now gives 

/""° dv 

m = vr - 2 / . , (8.37) 



^1-2^(6) 



where 



^^^^ - 2(^.-2)W2p2j,.-2 ' (8-38) 

and the turning point Uq is now described by the positive root of 

1-2^1^^-2-14^ = 0. (8.39) 

At infinite impact parameter, the scattering angle vanishes and so does ^{b). Thus for 
large impact parameters, the integral in Eq. (8.37) can be expanded in ^, giving 

eib) = 2 m I' du^^^ + oie) = ^"~'^r(!!/V'^^'^ ^^^^ + ^^^'^ • ^'-'^^ 

Making use of Eq. (8.38), we sec that at large impact parameters b, the momentum transfer 
(8.33) is given by 
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^^ = V|^^, (8.41) 



in which 



9^-^ = (- 



An J \ J Tl" ^ 



For finite impact parameters, no divergence appears if we simply set z/ = 3. In this case, 
the integral in Eq. (8.37) gives 9{h) = 2tan~^ ^{b), and with z/ = 3 in the definition (8.38) of 
C{b), one finds the well-known result for Rutherford scattering, 

i/ = 3 : 

Taking into account our results for at z/ = 3 and arbitrary but finite b, Eq. (8.43), and 
at large b but arbitrary u, Eq. (8.41), we arrive at an interpolation formula 

<lKb)=^p \,^_,\^.^.., (8-44) 

which is valid for u slightly above 3. Since [g^(6) — q]{b)] vanishes in the limit — > 3 for 
finite values of 6, while the integral of b"^"^ [q'^{b) — qj{b)] over very large b values vanishes, 
we conclude that, in the limit — > 3, we may replace the exact g^(6) by the interpolation 
function q]{b) in computing the momentum transfer integral (8.32), with 



1/^3+ . 

Vpb I da^, = — 1^.-2 4p' / , (8.45) 



in which 



dbb^- \_.,._: , ,_.._4 - (8-46) 



/o b^-^ + bt 

This integral is akin to the previous integral (7.21), and a similar evaluation gives^^ 



^^Namely, one sets 6^*^ ^ = ftg*^ ^ x to express / as a contour integral that gives the discontinuity 
of a;(3-^)/(2'^-4) . The contour may be then opened up to enclose only the simple pole at x = — 1, 
which gives the result (8.47). 
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I ^ iKl / = + 0(z/ - 3) . (8.47) 

Placing this result in Eq. (8.45) and using the definition (8.42) of 6o, we obtain, after a little 
algebra, the i/ — > 3 limit 



1 

— ^ -7 



.48) 



27r P \ 4p2 y 27r 

The pole a,t u — 3 that appears here refiects the long-distance, infra-red divergence that 
appears when u approaches 3 from above. Note that for simplicity, we have written the 
result for like charges, CpCf, > 0; otherwise this product should be replaced by \epeb\. 

2. Classical Coefficients 

We turn now to compute the transport coefficients when the classical cross section (8.48) 
is placed in the general result (8.30). To do this, it is convenient to use the velocity variables 
Vp and Vft, with p = mpj, Vpj, and Vpb — — v^. The classical cross section (8.48) produces 
the factor p-ip2(''-3)/(i^-2) 

vpb J da-,q'^vlp\ (8.49) 
As Eq. (8.30) shows, the coefficient involves ^pb/vph times this factor. Since 

and an integration by parts makes the derivative act upon the distribution function fb{Pb), 
we have, effectively, 

d 

>Pi,mb^b- (8.51) 

OVb 

Hence, for classical scattering, Eq. (8.30) may be expressed as 

To reduce this expression, we use the integral representation 

vlr^ = [^) I - exp |--/3,m.^J,5| . (8.53) 

With 
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(8.54) 



and with Vp^ = Vp — v^, we may interchange the integrals and complete the square in the 
Gaussian integral to get 



v—2 I 



z/ - 2 Vp ■ Vfe 



= rife 



2v-4 



V — A u 2 



1 



r 



exp{-iA-5-JYT7} {-^r^> l} • (8-55) 

We use this expression to evaluate the general formula (8.30) using the classical cross section 
(8.48). The result is simphfied by the variable change u — s s)"^ , and its nature clarified 
by introducing the squared Debye wave numbers for the particles of species 6, 



nl = (3b el Tib . 



Thus 



^ll> I _ ^ 



I/-3 



-7 



2 / 



1/2 r 



8m„h 



^ \2u-4) 



J^'diiii-V2^,(3-)/(-2) (1 -^,)K-3)/(2.-4)exp |-i/36m6^;>| {-^«, 1 



(8.57) 



To facihtate the comparison of this u > 3 result with the result (7.28) for i/ < 3, we first 
note that, in the limit in which u approaches three, 



1/ 



^ii;^) r(i) 



l + (^_3) 



with 



Using 



^[l-(ln4 + 7)(^-3)] + 0{u-3)\ 



i/^3: — ln4 = ^— 4^=^-'^) + C)(i/-3). 

u — 3 u — 3 



(8.58) 



(8.59) 



(1 _ ^)K-3)/(2.-4) ^ (1 _ ^)(.-3)/2 ^1 _ ^)(.-3)/{.-2) ^ 



(8.60) 
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and 



u-2 



i/-4 

we may now cast Eq. (8.57) in the form 



l + 2(i/-3) + 0(1/ -3)2, 



(8.61) 



A> 



r<ii> 



4^ 27r '^^ 
exp|-i 



-27 



(i-i) 



2mp6 1 ~ 'i*. 



{[l + 2(i/-3)] 1^, 1} . (8.62) 



IX. CLASSICAL RESULTS 

The > 3 contributions (8.62) we have just computed and the u < ?> contributions of 
Eq. (7.24) are each separately divergent in the i/ — > 3 spatial limit. However, it follows from 
the general principle of our dimensional continuation that their sum is well defined in the 
limit^°, combinations we denote by 

A' = lini{A",c + A^} , (9-1) 

and 

= lim{C^;> + CM. (9.2) 

It is worthwhile emphasizing again that this provides the leading and next-to-lcading 
contributions to the classical stopping power in the plasma coupling g. The lower-dimensional 
contributions (<) have several pieces, Eqs. (7.24)-(7.26). They consist of the sum of singular 
terms, Eq. (7.25), and regular pieces, Eq. (7.26). Note that it is redundant to employ a 



^'^In this regard, it is worth noting that as the spatial dimensions change, the physical dimensions 
of a charge change. This is made explicit if we replace . Now the wave number 

ji carries the dimensional change in the original definition of the squared charge. The fact that 
the singular pole terms cancel implies that the complete result is independent of what precise 
numerical value is taken for the arbitrary wave number ji; the result is insensitive to the way in 
which the squared charge is extrapolated away from three dimensions. This is the analog, within 
our dimensional continuation method, to the renormalization group invariance of quantum field 
theory. 
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classical subscript in the u < 3 case, as we did for the u > 3 piece, since all contributions in 
dimensions less than three are purely classical. 

Since the regular piece is finite and has already been reduced to its simplest form, let us 
concentrate on the sum of the singular pieces, and write 

= li^S { A".c + A",s} , C-g = Ihn {Ci% + Cl[<} . (9.3) 

The subscript S on the left-hand side of Eq. (9.3) should not be taken to indicate that the 
limit of the sum is singular, but only that the result comes from the well-defined sum of 
individually singular pieces. The — > 3 limit is well defined, as the respective pole terms 
cancel. The first term was derived in (8.62), and is dominant for i/ > 3, while the second 
term is dominant for < 3 and is given by (7.28), which we repeat here for convenience: 

Some computation shows that the u ^ 3 limit of the sum of Eq's. (8.62) and (9.4) gives 

1/2 



u 



(9.5) 



Since the classical energy loss functions A'^^ and C^'g contain the complete contributions 
for z/ > 3, the only additional finite part is that which comes from z/ < 3, and so the complete 
energy loss functions for the species 6 plasma particles in the classical case are given by 

= Al, + , = C^S + ^£ , (9-6) 

in which the two contributions to each function are given by the results (9.5) that we have 
just dealt with and the previous results (7.26), which we repeat here for convenience: 

An PbVp 27: J _i Ptotai(^'pC0s6') \ J { cos 



(9.7) 



Following the methods used to show that the result (3.2) is independent of K, it is straight- 
forward to verify that the sum (9.6) of Eq's. (9.5) and (9.7) is also independent of the 
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particular value of K. Moreover, the reflection symmetry F{—u) — F*{u) guarantees that 
the result (9.7) is real. 

The classical result apphes to the t'p — > hmit of the energy loss (large velocities are 
inconsistent with the classical limit). This limit of (9.5) entails elementary u integrals^^, and 
one finds that 



: 



In ( ^'^K^ ) + 27 



-^^m,vl^ + 0{vl). 



(9.8) 



This result, added to the small velocity limit of the regular part (7.34) derived in Sec. VII, 
produces 
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In ( A ''''''' 



Kr 



rrib 



IGtt nipb 



+ 2+27 



2 

15 



Pm^vl - I E 4 /^cme + ^ 4 iPcm^f'] ' I + 0{vl) . (9.9) 



Note that the arguments of the logarithms that appear here involve a small factor that 
is essentially the plasma coupling parameter g ~ (5e^Ko/^'K- Hence these logarithms are 
negative numbers that are large in magnitude. 
In a similar fashion, we compute 
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(9.10) 
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which, added to the small velocity hmit of the regular part (7.35), produces 



rrib 
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+ (9.11) 



Finally, we note that Bb — C"^ — Ab{'^/ l^bVp) has the leading small velocity hmit 



Vp 

B? 



47r 



1/2 

2T:(5bJ 3 



y/- 4 



In 



167r 



+ 2+27 



(9.12) 



X. QUANTUM CORRECTIONS 
A. Quantum Scattering in the Born Approximation 

No dimcnsionless parameter can be formed from the basic quantities 6^65 , nipb , and Vpb 
that describe the classical scattering. Hence dimensional analysis determines all of the 
classical result (8.48) except for the purely numerical factors that it contains. Quantum- 
mechanical scattering, on the other hand, is richer in that it involves the dimcnsionless 
parameter 

where Vpb = p/rUpb is the velocity of the projectile relative to a particle of species h in the 
plasma. We have previously made use of this parameter in the v — 2) limit. The factor 
(p/tt/i)'^'^"^^ that disappears when 1/ — ?> \s introduced to make ripb dimcnsionless when the 
spatial dimensionality u is extended away from v — 3.^^ The quantum-mechanical exten- 
sion of the classical squared momentum-transfer cross section (8.48) involves a dimcnsionless 



Apart from three dimensions, the Coulomb potential, being the i^-dimensional Fourier transform 
of behaves as l/r^^~'^\ Hence /r^^"'^^ has the dimensions of energy. Since Planck's con- 

stant h has dimensions of momentum times distance or, equivalently energy times time, an inverse 
distance has the dimensions p/h while hv has the dimensions of energy times distance. We conclude 
that Eq. (10.1) does indeed define a dimcnsionless parameter for arbitrary spatial dimensionality 
V. The additional factor of tt in {p/Trh)^'^~^^ is introduced for later convenience. 
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function of the dimensionless parameter rjpb. This parameter, which describes short-distance, 
quantum-mechanical effects, cannot appear to leading order in the long-distance physics 
which is involved in the u < 3 process evaluated in Sec. VII. It can, however, enter into and 
correct the short- distance scattering process with which wc arc now concerned. Although 
the quantum-mechanical description that we now turn to is strictly outside of the derivation 
of the Boltzmann equation from the classical BBGKY hierarchy, it is clear from physical 
grounds that quantum mechanics must be employed when ripb is small, which formally cor- 
responds to a limit in which Planck's constant h becomes large. In the limit in which the 
projectile velocity Vp becomes small and the rj parameter is large, the previous classical limit 
must be employed. Otherwise, a quantum-mechanical treatment of the scattering in v > 3 
must be made. 

We turn now to evaluate the scattering factor (8.15) multiplied by when rjpb is small 
and the quantum-mechanical Born approximation result is appropriate. This gives the ex- 
treme quantum-mechanical limit that applies for very high projectile velocities. We shall 
soon bridge the gap between the quantum Born and the classical results. In the Born ap- 
proximation, 

T=^, (10.2) 



and so 



We use 



J {2nhy 



= sin2^/2, (10.4) 

and express the momentum integration volume in (hyper-) spherical coordinates, with an 
implicit integration over all the angles save for the polar angle 6, to write 

d^p' = rupbp'^"'^^ d{p'^l2mpb) sin'^-^ 9 dO , (10.5) 

where ilj/_2 is the solid angle of a — 2 dimensional sphere. Setting 9 — 2x gives 



The integral which appears here has the value {v — 3)~^ + Oiy — 3), as one can show by 
dividing it into two parts with a suitable partial integration, or by expressing it in terms of 
the standard integral representation of the beta function. Hence 



Vpb j 



.B „2 _ (epCf,)^ TUpb ( V'' ^^^^ 1 



Again the pole which appears here reflects the infrared divergence when u approaches 3 from 
above. 
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B. Full Quantum Correction 



To fill in the region of arbitrary rjah values, we consider the weighting of the squared 
momentum transfer with the difference between the complete and first Born approximation 
cross sections, 



This integral of a cross section difference is well behaved in the limit ly ^ 3. The pole at 
V = 3 produced by the integral involving the full cross section is canceled by an identical pole 
in the integral involving the Born approximation. This is because these poles come from soft, 
infrared physics corresponding to long distances where the potential is weak. The divergence 
behavior leading to the poles is produced entirely by the first Born approximation term of the 
full cross section da which is then canceled by the subtraction of its Born approximation da^. 
However, we cannot simply set the spatial dimension z/ = 3 in the cross section difference 
that is the integrand in the integral (10.8) because in three dimensions the Born and full 
Coulomb cross section elements are identical, da = da^. The u ^ 3 limit of the integral 
is not the integral of the z/ ^ 3 limit of the integrand. The integral does not converge 
uniformly at small scattering angles, and the order of the limits cannot be interchanged. 
What we shall do is to implicitly assume that the spatial dimensionality is slightly greater 
than three so as to regulate the theory^^. Then we shall make a (implicitly generalized) 
partial wave expansion for the cross section difference (10.8). High partial waves with / ^ 1 
correspond to large impact parameter scattering where the effect of the potential is weak 
and the first Born approximation becomes exact'^^. Thus the subtraction of the first Born 
approximation within a partial wave decomposition yields a partial wave sum that converges 
at large / values and hence gives no pole at z/ = 3 when the physical limit of three dimensions 
is taken. Thus all we really need do is to express everything in terms of partial waves in 
three dimensions and subtract the Born approximation in the partial wave summand. In 
this way we may exploit some clever mathematics of Lindhard and Sorensen [29], but in a 
manner which justifies its use. It should be emphasized that the Born approximation must 
be subtracted in the partial wave summand before the sum is performed — the separate 
sums do not converge at large I. 



^^Other infrared regularizations may be used, such as the replacement of the Coulomb potential 
with a screened Debye potential, since only a potential logarithmic divergence is to be avoided. 

^^At large angular momentum, the effective centrifugal potential is much larger than the Coulomb 
potential, justifying the treatment of the Coulomb potential by the first Born approximation. 
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Although we always have in mind that the Born term is to be subtracted, for simplicity 
we shall omit this explicit subtraction in some intermediate steps. We use the standard 
partial wave decomposition of the scattering amplitude, 

da = dQ2\f{0)f , (10.9) 

with 

^ oo 
1=0 

The cross section weighted integral of — (1 — cos^) may now be evaluated using 

{21 + 1) cos9Pi{cos9) = (/ + 1) Pi+i(cos^) + i P,_i(cos^) , (10.11) 
and the orthogonality relation 

J {21 + 1) Pi, {cos e)Pi{cose) = at^s^, ■ (10.12) 

These give^^ 

/oo 
da g2 = 277^2 ^(Z + 1) {2 - e^^l-^'-^t'+i)! - e-^^l^'-^o+DJ } . (10.13) 
1=0 

For the Coulomb potential^^, 

r(Z + l-^r/) ' ^ ^ 

where the phase (j) is independent of and r/ is the generic quantum parameter. For the 
specific p-h system that we consider, r] —>■ rjpi,, where 

_ epei,mp, _ e^e, 



Appendix D explains how the cross section averaged momentum transfer (10.13) is simply related 
to the classical limit. 

^^This formula for the Coulomb partial wave phase shift is derived in many graduate level quantum 

mechanics texts. See, for example, Gottfried [30], Sec. 17, Landau and Lifshitz [31], Sec. 36, or 
Schwinger [32], Sec. 9.2. The connection of the Coulomb phase function to that of the Debye 
potential, and the recovery of the Coulomb result in the infinite screening radius limit, has been 
presented by Brown [33] using determintal methods and Jost functions. 
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Using r{z + 1) = zr{z), a little algebra, and subtracting the Born approximation, we find 
that [34] 



1=0 



+ 



l + l + ir]ph l + l-ir]ph ^ + 1 



(epebV 1 



47r V, 



— 2[Re^jj(l + ir]pb)+l] , 



(10.16) 



pb 



where t/j{z) is the logarithmic derivative of the gamma function, t/j{z) — r'{z)/r{z), and Re 
denotes the real part. Note that, by combining denominators, we may write 



RexPil +ir)) = y- ' 



-7 



;i0.17) 



To check the partial wave method that has been employed, we consider the limit in which 
Tjpb becomes large, in which case the full quantum cross section becomes equal to its classical 
limit. 



Vpb J {dapb - da^^) Vpb J {da'^^ - da^^ 



Indeed, using 



— > oo : 



'il){l + z)^\az + 0{z-^) 



;i0.18) 



(10.19) 



on the right-hand side of Eq. (10.16), we find that 



Vpb J {da^^ - da^^) q^ = 



{epCbf 1 



47r Vpb 



In 



( ^ 



eft 



yAnhvpb 



+ 27 



;i0.20) 



This is precisely the u ^ 3 limit of the difference of Eq's. (8.48) and (10.7), confirming the 
validity of our use of the partial wave expansion^^. 



37 Although we are not able to explicitly compute the cross section difference (10.16) for > 3, it 
easy to make a model of the mathematical expression for this case, a model that gives the essence 
of the ly > 3 behavior and which reduces to the correct ly ^ 3 limit. We set 

Vpb J {dapb - dapb) 1 = ^ [-^) ^ Finpb, u) . 

The overall factors here have the correct dimensionality of velocity times squared momentum times 
length to the power v — 1. We combine denominators as in Eq. (10.17) with a fractional power 
of 1/A; chosen, as we shall see, to give correspondence with previous results. Thus we choose the 
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It is convenient to refer the total cross section integral of the squared momentum transfer 
to the classical cross section, not the quantum Born approximation. Thus, we subtract 
Eq. (10.20) from Eq. (10.16) to obtain the purely quantum mechanical correction 

Vpb J {dapb - da^,) = {2 Re V^(l + tTjpb) - In 4b} ■ (10-21) 

Equation (8.30) gives the explicit form for the energy loss functions for the plasma species 
b in the u > 3 region that we are now considering. For convenience, we repeat this formula 
here: 

In view of this general formula and the expression (10.21) for the difference of the complete 
quantum cross section and its classical limit, we may write the v > Z result for the energy 
loss functions in the general case as 



model function 



n,,.)^-i:(i) 



oo 



fe=l 



^ k I ^ \ k J k"^ + ri^ 

k=i ^ ^ k=l ^ ' 

These formulae reduce to Eq. (10.16) in the limit u ^ 3. Of more interest is the character of the 
77 — oo limit. To obtain this limit for v slightly greater than 3, we note that the first sum in the 
second cquaUty above defines a zeta function, and with ({s) = {s — + 7 + 0{s — 1), we have, 
for 1/ — > 3, 

E/l\U-2j 1 
[k) =;7T3 + i+7- 

k=l ^ ' 

The second sum, in the large 77 hmit with v slightly above 3, may be replaced by the integral that 
provides the asymptotic value 



"-2 I fc^ 1 (^) 



,yU-27 + 1 . 



Hence, in the — 00 limit, 

^) = ^73^ {r?(^) - 1} - 7 , 

and we see that, using the definition (10.1) of our mathematical model reproduces, in this limit, 
the difference of the result (8.48) for the classical scattering integral and the result (10.7) for the 
Born approximation scattering integral (with the neglect of terms that vanish when z/ — > 3) . 
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Al = Al, + ^ , = eg + , (10.23) 

where ^g"^ and C^'^ are the z/ > 3 classical results given in Eq. (8.62), while A^'^ and c''^^^ 
are the quantum mechanical corrections to this classical result, the results given by inserting 
the correction (10.21) into the general formula (10.22). Explicitly, 

f_VVj^_jl (10.24) 

1^ PbrrihVpVpb ) 

In accordance with our general principle of dimensional continuation, to find the leading 
and next-to-leading contributions of the stopping power, we must take the i/ — > 3 limit of the 
sum of the v > 3 piece (as calculated above) and the u < 3 piece. The complete quantum 
result in three dimensions is therefore provided by 

Ab = lim {Al +At}^A^b+ , (10-25) 

and 

Ci^ = lim {Cl^> + Ci^< } = Cj}"" + C"^^ , (10.26) 

I/— >3 

which, with the aid of Eq.'s (10.23) and (9.1) - (9.2), we have written in the form of a purely 
classical piece plus a quantum correction. We have already calculated the classical term 
(9.6) and found that the potential divergences cancel, so we turn to simplifying the quantum 
piece (10.24). Namely, the integration variable can be changed to the relative velocity, and 
the angular integrations performed. This gives 

expj -PbTUbivp '"p^^^\^i^^^^y^y^^{^ (^bmbVpVpb) 

+ exp I -(3bmb {vp + Vpb) | | ^^^^^^^^^ + p^^mbVpVpb) ' } ' 

(10.27) 

This expression provides a small correction when the integration is dominated by regions in 
which the quantum Coulomb parameter rjpb is large so that the scattering is nearly classical. 
When the effective rjpb values are of order unity, then a detailed evaluation of Eq. (10.27) is 
called for. But there are some limits in which Eq. (10.27) simplifies. 
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C. Simplifications and Asymptotic Limits 



One simplification appears for cold plasmas, that is, plasmas for which the thermal speed 
Vt = \/3/(3bmb can be neglected. This is described by the formal limit (3b ^ oo. In this limit 
the first exponential in Eq. (10.27) sets Vpb = Vp in the factors that multiply it, the second 
exponential gives a negligible contribution, and so Eq. (10.27) becomes. 



Vb 



oo 



47r 2(3bVp 

exp{-iw}{^.l 



47r 2(3bVp 



l^{2R.,(l+.,,)-ln,^}{^.l}. 



(10.28) 



where 



€p6b 

47rhvn 



(10.29) 



The other simplification appears when the thermal velocity or Vp or both are large in 
comparison with ebCp/ATrh. In these cases we may use the small rjpb limit which Eq. (10.17) 
reveals to be^^ 
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Reip (1 + IT]) = -7 + 0(?7^) 
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It is difficult to implement this limit directly in Eq. (10.27). It is much easier to return to 
the starting point (10.24) and evaluate it in the manner of the evaluation of the classical 
energy loss functions. Since = m^^Vp^, the starting point Eq. (10.24) involves 



{2Re^l:{l + ir]pb)-\nr]l^] 
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^^Because of the small Vpb integration region in Eq. (10.27), the formal order r/^ error in this limit 
is actually of order 77^ in 
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We first express 
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lim - 

6-»o e 



-pbmbVpf, 1 - 1 



;i0.32) 



and then exponentiate tlie terms involving Vpb using 



-I3bmbvl^ 



I f^ds ^ r 1^ 2 



with p = 3/2 - e, p = 3/2, and p = 1/2 - e, p = 1/2. Writing Eq. (8.58) in the form 

_i_ = _L[i_(log4 + ,)e] 

and using zV{z) = r{z + 1) to also obtain 

1 2 



;i0.33) 



10.34) 



r(3/2 - e) y/^ 



l + (2-log4-7)e] + 0{e^) , 



(10.35) 



we find that 



1 



{2 Re-i/^ (1 + irjpb) - InT^JJ {l , /^bmbvJJ 

\ 3/2 



= - ( 2^^"^M y ciss^/^exp |--/?6m6'ypj 
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Placing this representation in Eq. (10.24), interchanging integrals, performing the resulting 
Gaussian integration, and making the variable change previously used, 

u 



1-u 



(10.37) 



yields 
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c2^2 / \ 1/2 ^1 



J duu ^^'^ exp I — -/SfemfeVpti 
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The rjpb limit is formally the large h limit. This is the limit in which quantum 
uncertainty rather than a classical turning point sets a minimum distance scale. This may 
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be brought out explicitly if we replace the previous combination (9.3) of potentially singular 
parts by 



b.S 



•^h 1 



^llQ _ ^Uc 1 /^llAQ 



'b,S 



;i0.39) 



The complete functions defined by Eq's. (10.25) and (10.26) now read, in view of Eq. (9.6), 



where the regular parts coming from the u < 3 contribution were defined in Eq. (9.7). Adding 
Eq's. (9.5) and (10.38) gives the quantum regime limit 
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2 \mpb 2mpb 1 — u J Jl J I J 



(10.41) 



To simply compare quantum formula (10.41) with the classical formula (9.5), we neglect 
mass ratios. Then we see that the classical cutoff length /S^ (epeb/47r) in Eq. (9.5) is replaced 
by the quantum length h \J j3h/ 2mpb here in Eq. (10.41). 

Equation (10.41) is the high velocity limit in the quantum-mechanical sense that 
\Vpb\ ^ 1) But this limit entails no restriction on the kinetic energy mhVp/2 relative to the 
temperature = l/Pb- When the kinetic energy is also large in comparison with the 
temperature, a limit that we shall simply denote as t'p — > oo, the exponential damping in 
the integrand of the integral (10.41), which emphasizes the u — region, allows us to set 
1 — w — > 1 in the logarithm and extend the upper integration limit to u = oo. In this way, 
we obtain 
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(10.42) 



'''/Sb^^p' ^"'^ J 47r (5bVp 

This is to be combined with the large velocity hmits (7.41) and (7.47) of Ab^^ and C"^ . The 
coefficients Ab in Eq's. (10.42) and (7.41) are both dominated by the electron contribution, 
a contribution mb/me larger than that of an ion h. To bring out the nature of the result, it 
is convenient to use the squared electron plasma frequency uj^ — Kg/zSgme, and we have 



oo 



(10.43) 
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On the other hand, all plasma species contribute to 



Vp ^ oo 



C^' = ^ In f ^-^] . (10.44) 

Since Ah{'^/ Ph'^p) behaves as l/v^, it vanishes more rapidly than C" for large Vp [of relative 
order 1/ ISbTUpVp] and so the general connection Bb — C" — Ab{'^/ l3bVp) gives 



— >• oo 



CpHh 



47r /^bVp V Hkb ) ' (10.45) 

The limit that we have just described applies to the situation in which the projectile 
velocity Vp is so big that irieVp » Tg. For the case where the projectile is an ion, this 
limit implies that the projectile kinetic energy Ep — mpVp/2 is even greater than a typical 
plasma temperature T by the additional large factor of mp/nie. The hmit in which Ep^T 
but yet rUeVp is not large in comparison with T is also of interest. In this case, so long as 
{Ep/T)^ [me/mp) > 1, Ae still dominates over the ionic contributions to Ah, but a detailed 
evaluation of Ae is required. In this intermediate case, the limit (10.45) still holds for the 
contribution of the ions in the plasma, while the contribution of the electrons in the plasma 
is of relative order {nieVp/Ty/'^ and thus may be neglected. 

XI. TRANSPORT EQUATION VALIDITY DETAILS 

Here we provide the detailed computation of the result (6.3) used in Sec. VI for the 
error of the Fokker-Planck equation as measured by different evaluations of the increase in 
transverse energy. Again we note that since the u < 3 contribution to the Fokker-Planck 
coefficient C^' is the same Lenard-Balescu equation that is used to evaluate the transverse 
energy in this region, the difference defining A;, is given by just the u > 3 parts. 



dEl, 



dt 



dElb 



exact 



(11.1) 



F-P 



with both terms computed from the scattering cross section formula that is equivalent to 
the Boltzmann equation as is described in Sec. VIII. This is the computation to which we 
now turn. 

The exact transverse energy change in a scattering with the initial and final projectile 
momenta Pp and p^, is given by 
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Ppl. 



Pp Pp - (Pp • Pp) 



;ii.2) 



2mp 2mpp2 

where we now append a subscript and write the projectile mass as rUp to avoid possible 
confusion. Here Pp = q + Pp, where q is the Galilean invariant momentum transfer, and so 

1 



[W ^ - (Pp • q)^] ■ 



(11.3) 



Since the relative velocity Vp^ is the only vector available to describe the Galilean invariant 
cross section, the tensor 

can only involve the tensors v^^v^^ and d'^K Prom energy conservation in the center-of-mass 
system. 



1 1 
Vp6 • q = 7^{^'pb + Vpfe) • mpb{w'pb - Vpb) - -(Vp6 - Vpb) ■ mpb{Vpb - Vpb) 
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= 
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(vp6 - vjfe) - - rupb (v^fe - VpbY 



q^ 



2mpb ' 



where 



1 ^ 1 



(11.4) 



(11.5) 



TTipb rrip rrib 

defines the reduced mass rUpb- Thus, by contracting with and with v^^v^^, it is easy to 
verify that, weighted by the cross section, we have, effectively. 
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)m2 ^2 
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pb pb 



Accordingly, we may write the transverse energy change (11-3) for scattering off particles of 
plasma species b as, effectively. 



AE 



q^ 



±b 



2mppl 



Pj - (Pp - (Pp • Vpb)^) 



(11.7) 



This appears in the general formula (8.9) which entails an integral involving 
J d^PbfbiPb) ~ J d'^^pb exp |-^/3fomfe (vj, - Vpbf^ 

PIT 

~ / sin^~^ 9 dO exp { I3bmb VpVpb cos 9 } , 
Jo 



(11.8) 
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with 



;ii.9) 



and the remainder of the integrand independent of the angle 9 between Vpb and Vp. Hence 
we encounter 



/3j,m(, VpVpt, cos d 



/ d9 sin^ e sin"-'^ 9 e^""*" ^ = - / sin^"^ ^ d 

Jo Jo Pbrnb VpVpb 

Jo Pbrrib VpVpb 



^ ^PbmbVpVpi, cos 6* 



with the last line following by partial integration. We thus have, effectively, 

{v — 1) cos^ 
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RecaUing the relationship (4.30) of the Fokker-Planck approximation to the rate of trans- 
verse energy increase and the connection (4.24) amongst the Ab, Bb, and coefficients, we 
have 



dt 



F-P 
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PbVp 



(11.14) 



Making use of the formula (8.30) for the Cjj^ and A^ coefficients, and referring to the general 
formula (8.9), we see that the previously defined error measure is given by 



■ fbiPb) Vpb / dapb < AE± - ^ p p 



(11.15) 



2mp PbTrib VpVpb '^rUp 

All the potential infrared singular terms, the terms involving a single power of q^, cancel, 
and there remains, in the limit v ^ 3 which now may be taken, 

d^Pb 



fb{Pb) 



Snip m^f^Vpb 
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Vp • Vp6 
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J dapb{q^f. (11.16) 
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The cross section weighted integral of {q^Y that appears here may be evaluated, for example, 
by inserting an extra factor of — Ap^ sin^ 9/2 in Eq. (10.6) restricted to u — 3. The result 
is that 



I dapb {qr - ^ 



m: 



pb ■ 



We write 



^Ib 



d 1 
dvb Vpb 
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(11.18) 



and integrate the velocity derivative by parts so that it acts on the distribution function 
MPb), giving, effectively. 
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dvh 
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Hence, 
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We write the factor 1/vpb in terms of a Gaussian integral, interchange integrals, complete 
the square, and change variables as in the computation of Eq. (8.57). This gives 



Att 2mr, 



which is the result previously quoted in Eq. (6.3) in Sec. VI. 



XII. RATE AT WHICH DIFFERENT SPECIES COME INTO EQUILIBRIUM 

Plasmas may be created that contain different species which are at different temperatures. 
This happens, for example, when a plasma experiences a laser pulse which preferentially heats 
the light electrons that have the larger scattering cross section. Here we use the methods 
that we have developed to compute the rate at which the various plasma species come into 
thermal equilibrium. 



A. Introduction and Summary 



We shall assume that the particles of two species a and h in the plasma arc individually in 
thermal equilibrium, but at different temperatures Ta and T^. We shall compute the leading 
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and subleading orders, as we have done throughout, of the rate d£ab/dt at which the energy 
density £a of species a changes because of its interaction with species b. Since 

where /a(Pa) is a spatially homogeneous thermal distribution at temperature = the 
Fokker-Planck equation (4.1) gives 

^ = (T„-T,). (12.2) 

in which 

Cab^ j ^^fa{Va)l3aVaA,{Pa). (12-3) 

Since the energy loss of one plasma species is another's gain, the rate of energy density 
transfer is skew-symmetric, 

dSab d£})a ^^n A\ 

~df~~~df' ^ > 

whence 

Cab = Cba (12.5) 

are symmetric coefficients. 

A plasma consists of light electrons of mass rUe = m and heavy ions of mass ^ m 
which we shall generically denote by M. Before plunging into the details of our computations, 
we review the well-known justification for assuming that the electrons and ions in the plasma 
are themselves in internal thermal equilibrium at the separate temperatures and Tj. To 
do this, we note that, as shown for example in Eq. (12.25) in the results below, the mass 
dependence of the rate appears predominately in 



Jmanih 

= -7 .3/2 ■■■■ (12-6) 

[niaTb + rribTa) 

Let us now use this result to compute the rate at which the electrons come into equilibrium 
with themselves. For this purpose, we imagine the very simple situation in which the elec- 
trons are in two pieces, one with temperature T^, the other with temperature T^. Then the 
rate at which these two pieces come into equilibrium is controlled by the factor 

1 



(7; + 76)^/2 mV2 



:i2.7) 
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If a similar partition of the ions into two parts at different temperatures were made, the 
parts would come into equilibrium at a rate controlled by the factor 

^ fl2.8) 



(r„ + n)3/2MV2 



Thus the rate at which the ions come into equilibrium amongst themselves is a factor ^Jm/M 
slower than the corresponding rate for the electrons. Now going back to our original problem 
of electrons and ions at different temperatures Tg and Tj, we see that if the temperatures 
are not greatly different, the rate at which the two species come into thermal equilibrium is 
controlled by the factor 

which is a factor oim/M smaller than the rate at which the electrons come into equilibrium 
amongst themselves and a factor y^m/M smaller than the equilibrium rate for the ions alone. 
Thus our work which treats the electrons and ions as being in separate thermal equilibrium 
but with two different temperatures and concentrating on computing the rate at which the 
ions and electrons come into thermal equilibrium is justified by the very large ion - electron 
mass ratio. 

Let us consider the case in which the electrons have come to temperature Tg and all ions 
species have equilibrated to a common ion temperature Tj. With electron and ion specific 
heats per unit volume Cg and Cj defined by dS^ — CgdTg and d£i — J^i^i — Ci^Tj, we define 
the the rate F at which the ionic and electronic temperatures come into equilibrium by 

j^{Te-Ti)^-r {Te-Ti) , (12.10) 

with 

r = Cei (- + -] ■ (12.11) 

The sum of Eq's. (12.44) and (12.57) [evaluated for K — Kg as Eq. (12.44) requires] give the 
limit 

Time < Teirii : 

where 

(12.13) 
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is the sum over all the squared ionic plasma frequencies. The overall coefficient of the 
Coulomb logarithm in the energy transfer rate (12.12) was obtained long ago by Spitzer [35] 
as described in his book [36]. However, our determination, as always, gives not only this 
coefficient, but also a precise definition of the value of the Coulomb logarithm, the constant 
under the Coulomb logarithm. 

The development in the subsequent sections follows the order that we have previously 
used in the stopping power work. 



B. Classical Results 



The decomposition (9.6) previously given for the classical contributions to the "S" and 
"R" contributions to the Ab functions gives a corresponding division of the temperature 
equilibrium coefficients. 



/7C _ /7C I ^< 

The result (9.5) for placed in Eq. (12.3) gives 



(12.14) 



^ab,S 



Paelnl f iSbrrihY^^ f d^Pa . / x 2 /"^ 1/2 / 2 1 



In 



4:71 rriab l-uj 



(12.15) 



The particle number density may be expressed 



as 



/(p) = n ( ^ ) ciV exp j -- ^mv" 



{27rh) 



2% J 



(12.16) 



The resulting Gaussian integration yields 



C-3--4^ /A^y/^ f\^^y, H:l:nu,r' 

An \2n J J, (P^rua + PbTUbuf 

In ( ^ K ^^^^ — 



47r ruab l-uj 



(12.17) 



To place this in a form that exhibits the symmetry under the interchange of the a, h labels, 
we change integration variables to 



s — 



f^brribU 2 
1 — u 



ab 1 



(12.18) 



where 
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2 _ Pama + Pbmb 



1 1 

+ 



+ 



(^arriaPbmb Pairia Pbirib nia rrib 
is an average squared thermal velocity. This gives 

hIkI 3 {(3amal3bmbf''^ 



;i2.19) 



47r ^ ^j3^rna + ^brrib 
In 



3/2 



T^^^TT^ 1+27-2 



;i2.20) 



The integrals that appear here provide a standard representation of the beta function, 



B{x,y) 



n^ + y) " Jo + 



ds ■ 



This is so because 



Jo ^ ^ r(5/2) 3' 



and, using 



we also have 



In s = lim 



POD 

Jo 



"^^riB/ll"^ [^(3/2) -^(1)] = ^ [2 -2 In 2], 



(12.21) 



(12.22) 



(12.23) 



;i2.24) 



and thus the evaluation 



1 Y^' (/5„m„/5,m,)V2 



In 



CaCb K 



47r AmabV^^ 



+ 27 



.27ry (/3am„ + /3bmb)'/' 
The long-distance, plasma screening correction A^^ presented in Eq. (9.7) gives 



(12.25) 



An ' 



ab,R ~ A _ f^a 



{2TThy 



+1 



fa{Pa) 

Pb{Va COS 9) 



— / dcos9 

27rJ_i Ptotal(^^aC0S6') 



F{va COS 9) In 



F{Va COS I 



Va COS 9 , 



where, we recall, the function F{u) may be expressed in the dispersion form 

/''LoLal(t') 



F(u) = - 



dv ■ 

u-v + irj 



(12.26) 



;i2.27) 



with the limit — >• 0"*" understood. The spectral weight is defined by 
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ptotai(^^) = y; Pc{v) , (12.28) 



where 



Pc{v) = nlv\l^e^v[-\Pom.v' \ ■ (12.29) 



We insert 



/ + 00 

dv5{v-VaCos9) (12.30) 
-oo 

in the integrand of Eq. (12.26), use the form (12.16) of the particle number density, inter- 
change integrals, and perform all the integrals save that involved in the insertion (12.30). 
Thus 

"''''^ 27r V 27r V 27r y/ 

dvv^ e^v[-\[Pama + P,m,]v'\ . (12.31) 

J-oo 12 J 27r Ptotal(^^) \ J 

Note that this formula exhibits explicitly the symmetry C^^^^ — C^an- 
Since 

F{v) - F{-v) - F{v) - F{vy = 2mptot.i{v) , (12.32) 
the In K dependence of C^^ ^ appears in 

= ^ j J dvv exp |-- + Am,] . j Ini^ + • • • 

= .y, ^^""^"^""^^^'l f 1 In + . . . . (12.33) 



(/3„m„ + /35m5)'/' V27ry 

Hence the sum (12.14) of Eq's. (12.25) and (12.31) is independent of the particular value of 
the arbitrary wave number K as it must be. 

Although a numerical computation is needed for the general evaluation of C^^, ^, it can be 
found when /3e?7?.e <S Am^, or 

TiiTie < Term , (12.34) 

where is a typical ion mass and Tj is the common ion temperature. This is the case that 
is usually of interest in applications. Since me/rrii < 10~'^, this constraint holds unless the 
ion temperatures arc very much larger than the temperature of the electrons in the plasma. 
As a first step, we write Eq. (12.31) as 
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r'"dvvexp\--P,meV^\-^^ — F{v)ln(^] . (12.35) 
This limit under consideration is formally equivalent to the limit mg — > 0, and so 

TiUle < TeTTli : 

E.Pi(^^)=Ptotal(^^)- (12.36) 



Hence, defining 

we have 

TiTTle -C Temj : 



Q,K = E-^e1,a, (12.37) 



— oo 



The integrand that appears here is analytic in the upper-half complex v plane. Remembering 
that we are to take the formal limit me — > first, for large v we have 

,2 



where 



2 



c^N ^ = ^ , (12.40) 
are the squared ionic plasma frequencies. If we take 

K' = 4, (12.41) 
then In {F{v)/K'^) vanishes for large v, and 

|v| ^ oo : 

m^(m)^-yA- (12.42) 



Thus at large \v\ the integral behaves as j dv/v . We add and subtract the corresponding 
contour integral over a semi-circle C at infinity in the upper half complex v plane. The 
original integral with the ir] prescription is equivalent to one over a straight line just above 
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the real axis. Adding the integral over the infinite semi-circle gives a closed contour integral 
enclosing no singularities which thus vanishes. There remains 

' [± = ± f\de = -l (12.43) 



and therefore the evaluation 

TiUle < TeUli : 



Qk = --^ (^X'W^l (12.44) 



2 27r V 27r ^ 

It should be emphasized that this evaluation is valid only for the choice K — Kg- 

C. Quantum Correction 

The complete coefficient is the sum of the previous classical result and a quantum cor- 
rection 

Ca6 = C + C^, (12.45) 

corresponding to the separation given in Eq. (10.25). Using the result (10.27) for Af^, we 
have 

(12.46) 

The relative velocity is defined by = — v^, and 

^a^ = 7^- (12-47) 

ATThVab 

The function ip{z) is the logarithmic derivative of the Gamma function T[z). We use the 
expression (12.16) for the particle number densities /a(Pa) and fb{Pb)- Making use of the 
relative velocity, we may write the resulting product of exponentials as 

exp /3„m„w^| exp /^brnftW^j = exp ^ 



exp <( iPafria + Pbrrib) 



PbTUb 



;i2.48) 



where we have made use of the effective thermal velocity defined in Eq. (12.19). We change 
the integration variables from , p;, to , Vab and perform the resulting Gaussian integral 
in Vq to obtain 
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r°° dii^ ( 1? 1 

I ^cxp|-^| {2Re^(l + .^.,)-lnr/^J . (12.49) 

We have written Eq. (12.49) in a form that makes its dimensions obvious: [k^ v\ = cm~^ sec~^. 

The hmiting behaviors of the result (12.49) are exhibited if we make the variable change 
v^jj — V2,C) which expresses 

.AQ_ 1 2 2 {f3amaP,m,y/' f 1 



2 (/3,m„ + /3„m5)3/2 V27r; 

^ 2°°^^Cexp|-^c} {2ReV^(l + %,C-'/') -ln77^,r'} , (12.50) 

in which 

Here the exponential damps large C, values, and so as far as evaluating limits are concerned, 
we can consider C, to be of order unity in the curly braces in the integrand in Eq. (12.50). 

The low temperature hmit corresponds to the low velocity hmit Vgh 0. This corre- 
sponds to the limit f\ah oo, which is the formal limit /i — > 0. This is the formal classical 
limit. Indeed, since 

X — > oo : 

{2Re^/;(l + ix) -Ina;^} ^0, (12.52) 

the quantum correction vanishes in the classical limit as it should. The major case of 
interest is the rate at which ions and electrons in a plasma come into thermal equilibrium. 
Because of the very large ion/electron mass ratio, mj/me S> 1, to a very good approximation 
1/J = Te/mg. The condition that f\e.i » 1 is thus equivalent to Te <^ (cjCe)^ me/(47r/i)^, or 
that the temperature is much less than the binding energy of a hydrogen-like atom. At 
such low temperatures, our assumption that we are dealing with a fully ionized plasma is 
generally invalid. Hence this classical limit is only of limited physical interest. 

The more relevant high temperature limit corresponds to the high velocity limit in which 
f\ah becomes small. Since ■0(1) = —7, we have 

?7a& < 1 : 

^AQ^_1 2 2 {Pafng pbmbY/^ fl_y^^ 
2''"''^(/3„m„ + /3,m,)3/2 [27rJ 



2 A I 2J t V C 



.^2 ..2 PamaPbrribY/^ f 1 1 , , fVabXX (.^ 
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The first correction to this result is of relative order 77^^. The case of interest is the 
electron- ion energy exchange rate where, using the notation Cj = Zje, 



2 e-me 



^Z2 27eV 



(12.54) 



is very small. Hence for the case of interest, the limit (12.53) suffices. This limit combines 
with the previous calculation (12.25) of C^^g to give, with sufficient accuracy. 



2 2 



1/2 



The complete coefficient is given by 



In 



hK 



;i2.55) 



C< -I- 
^ab,B. ^ 



^ab,S ' ^ab 



12.56) 



For the case of ion-electron relaxation, since me/ mi <^ 1, the electron mass can be 
neglected relative to that of the ion. Assuming that the ion temperature is not more than 
an order of magnitude larger than the electron temperature, TirUe <^ TerUi. With this 
restriction, the expression (12.55) simplifies, and the result may be expressed as 



2 Peme 1 

27r '^^ V 27r 2 



In 



-7 



(12.57) 
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APPENDIX A: SIMPLE EXAMPLE ILLUSTRATING THE METHOD 

Since the method used in this work is a novel one, we include in this Appendix a pedagog- 
ical, simple mathematical example that illustrates the basic idea. This is the computation 
of the behavior of the modified Hankel function Ki,[z) in the small argument z limit with 
the index u also small. The argument z will play the role of the small parameter in our 
work; the index u will play the role of the dimensionality except that in this simple Bcsscl 
function example we shall examine the region where u is near zero, not three. This example 
already appears in the preliminary account [9] of the new use of dimensional continuation, 
but it worth repeating here so as to have a clear, self-contained presentation. 

The Hankel function has the integral representation 

Although k is simply a dummy integration variable, it is convenient to think of it as a wave 
number or momentum variable. When z is small, exp | — | (^ + |)} may be replaced by 1 
except when one or the other of the factors exp{— 2; k/2} or exp{— 2;/(2 k)} is needed to make 
the k integration converge in the neighborhood of one of its end points. When u is shghtly 
less than zero, the integral (Al) is dominated by the small k, "infrared or long-distance" , 
region. In this case, only the exp{—z/{2 k)} factor is needed to provide convergence, and we 
have 

u <0 : 




The variable change k = z/{2t) places this integral in the form of the standard representation 
of the gamma function, and we thus find that the leading term for small z in the region i/ < 
is given by 

i/ < : 

where 7 = 0.5772 ■ ■ ■ is Euler's constant. Note that the second line describes the behavior 
for z/ < near z/ = including the correct finite constant as well as the singular pole term. 

When V is slightly greater than zero, the integral (Al) is dominated by the large k, 
"ultraviolet or short-distance" regions. In this case, only the eiq){—zk/2^ factor is needed 
to provide convergence, and we have 
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1/ > : 

The integral again defines a gamma function, and so 
1/ > : 

K{z) ^ ^ {^y\l - u^) , (A5) 

with again the rcsuh containing the correct finite constant as weU as the singular pole term. 

The result (A3) for u < can be analytically continued into the region z/ > 0. In this 
region it involves a higher power of z than that which appears in the other evaluation (A5), 
and hence this analytic continuation of the leading result for i> < into the region u > 
becomes subleading here. Similarly, the result (A5) for i/ > may be analytically continued 
into the region ly < where it now becomes subleading. An examination of the defining 
integral representation (Al) shows that these subleading analytic continuation terms are, in 
fact, the dominant, first-subleading terms. For u > one term is leading and the other 
subleading, while ior u < their roles are interchanged. Thus their sum 

KA^) ^ ^ { (i) [1 - H - [1 + H } (A6) 

contains both the leading and the first subleading terms for both u > and u < 0. In 
the fimit i/ — > the ( "infrared" and "ultraviolet" ) pole terms in this sum cancel, with the 
variation of the residues of the poles producing a logarithm, yielding the familiar small z 
result 

Ko{z) = -ln{z/2)-j. (A7) 



^^For example, subtracting the leading term (A3) for v < from the integral representation (Al) 
gives 



-zk/2 



1 



-z/{2k) 



For > I' > — 1, the integral on the right-hand-side of the equation converges when the final 
exponential factor in the integrand is replaced by unity, the z ^ limit of this factor. Hence this 
final factor may be omitted in the evaluation of the first subleading term. A partial integration 
presents the result as 

Jo 

whose evaluation gives precisely the analytic continuation of the leading term (A5) for u > 0. 
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It must be emphasized that the correct constant terms [In 2 — 7] are obtained by this method 
in addition to the logarithm — In 2; which is large for small z. The result (A6) was derived 
from the analytic continuation of results that were easy to compute in one or another region 
where either "infrared" or "ultraviolet" terms dominated. This is the essence of our method. 
Of course, the general result (A6) could be obtained by a more careful computation of both 
the leading and first-subleading terms in either of the separate > or < regions 
as was done in the previous footnote. In the work of the present paper, however, such 
an extraction of the subdominant terms would be very difficult indeed, although perhaps 
possible in principle. 

APPENDIX B: CONVERGENT KINETIC EQUATIONS 

Convergent Boltzmann transport equations have been discussed by Frieman and Book 
[11], Weinstock [12], and by Gould and DeWitt [8]. These are equations of the usual Boltz- 
mann equation form, but with modified collision terms on the right-hand side that account 
for both the hard, short-distance collisions and the soft, infrared, long-distance scattering 
that is modified by the background plasma medium. Just as in our work, these equations are 
valid only to leading order in the plasma density. In this appendix, we shall describe these 
convergent kinetic equations and then sketch how they are equivalent to our method which 
uses dimensional continuation. But before passing to these details, we should again note 
that our method gives only the correct leading order terms with no spurious higher-order 
terms. This simplicity of computation is to be contrasted with the kinetic equation approach 
which does yield spurious higher-order terms that must be identified and discarded to obtain 
a consistent, leading-order result. 

Let us first recall that the collision integral in the Boltzmann equation is of the generic 
form (8.2), which may be written as 

The papers cited in the preceding paragraph work in i/ = 3 spatial dimensions and write the 
total coUision term as 

CT'''\Pa) - C^-'iPa) + CTiPa) , (B2) 

where each of tlic two collision terms on the right-hand side have the generic form given in 
Eq. (Bl). 
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The first part C^b^^{Pa) accounts for Coulomb scattering taken to all orders with the first 
Born approximation subtracted so as to avoid double counting since it is contained in the 
second C^^{pa) term. As Gould and DeWitt note, and as we have spelled out in some detail 
in the discussion of the hard scattering corrections to the cross section weighted momentum 
transfer integral (10.8), the treatment of this hard collision contribution requires some care 
since, in = 3, the Born and exact Coulomb scattering cross section elements are identical. 
Gould and DeWitt regulate this contribution by taking it to be the scattering for a Debye 
screened Coulomb potential. They write 

(B3) 



I y hard I *^ | 1 2 



7^(1) 

-L D 



where is the full, all-orders amphtude for the scattering on a Debye screened Coulomb 
potential, and T^^^ is the first Born approximation to this amplitude. The squares of the 
corresponding amplitudes are subtracted to avoid the double counting mentioned above when 
the second collision term C^^{pa) is included. The amplitude squared iTop is asymptotic 
to the exact Coulomb scattering amplitude squared at large momentum transfer 

q = Pa - Pa = P6 - Pfe ' (B4) 

but its behavior for small does not include the correct soft physics which entails frequency- 
dependent, dynamical screening. The Debye screening makes the separate contributions 
of each of the two terms in |T'ii3-''dj2 ^j^g energy loss finite in the infrared region. The 
subtraction of IT^^"*!, however, makes the contribution of the difference (B3) to the energy 
loss finite when the Debye screening is removed, when Ko — ^ 0. But the price paid for this 
is a spurious unwanted contribution of — |Tn^-*p in the ultraviolet. This piece will take care 
of itself upon adding the correct infrared physics provided by |T^°^*p, where T®"^* is given by 
the first Born approximation to the dynamically screen Coulomb amplitude 



nsoft ^a^bfl' 



q^e{q/h,AE/h) ' 
with the energy change of the scattering process being 



(B5) 



A£;=P^--P^ = -P^-P^. (B6) 
2ma ^rua ^rub 2mb 

The dielectric function e(k, c<j) is that given by the random phase or one-loop, single ring 
approximation^*^ {7.2). Note that — |T^^'*p vanishes in the ultraviolet since T^°^^ 



^''Xhe soft contribution C^^^'(Pa) is also discussed in Section 46 in the Physical Kinetics volume 
of the Landau-Lifshitz series [21]. 
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and To ' are asymptotic as q ^ oo. Therefore, the addition of to |T''i^'ip cancels 

the unwanted ultraviolet part of mentioned above, leaving only the correct large 

behavior of |TdP- To reiterate. 



I ^^converge i 



hard 1 2 



IT 



(1) 



soft 1 2 



+ \T 



soft 1 2 



has both the correct ultraviolet and infrared behavior. The subtraction of IT, 



(1)|2 



(B7) 



docs not 



just avoid double counting. It also, on the one hand, removes the arbitrary dependence 
produced by |TdP, and on the other hand, removes the large contribution of jT*^"^*!^. 

This convergent kinetic theory approach is certainly valid, but it entails spurious higher- 
order corrections in the plasma density that must be discarded after calculations have been 
performed as Gould and DcWitt correctly do. In contrast to regularization (B3), however, it 
is simpler (but ultimately equivalent) to regulate the hard scattering contribution by contin- 
uing it to a spatial dimensionality u that is slightly above z/ = 3 and using the dimensionally 
continued exact pure Coulomb scattering amplitudes. With this regularization 



hard I ^ 



rp(u>3) 



.(r/>3)(l) 



(B8) 



and the phase- space collision integrals are also extended to > 3. This regularization 
automatically entails no additional, spurious, higher-order terms. 

To establish the connection of the convergent kinetic approach with our method using 
dimensional continuation, we add and subtract the first Born approximation to the scattering 
with a Debye screened Coulomb potential so that the soft collision term appears as 



CT{Va)=CT{Va)+C'^iM). 



~(soft 



(1) 



(B9) 



The first term here is given by the Boltzmann collision integral with the squared scattering 
amplitude replaced by 



\T 



soft 1 2 



(gV^^) e(q//i, A£;/n) 
where the second term is produced by 



(BIO) 



(1) 



1 



(Bll) 



Since It*^"^*] vanishes rapidly for large momentum transfer q, the collision integral C^^^^^ may 
be replaced by the Lcnard-Balcscu form which is employed in our method. This is explained 
in detail in the following Appendix C. 
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We now give a quick proof that the 'convergent kinetic equation' is equivalent to our 
method of dimensional continuation provided that the kinetic equation is solved in a consis- 
tent fashion and spurious, higher-order terms are discarded. To do this, we take the hard 
Coulomb scattering part in Eq. (B2) to be defined by the dimensional continuation with 
u > 3, as indicated in the discussion of Eq. (B8). We take the soft part in Eq. (B2) to be 
divided as in Eq. (B9) with the first term with the over bar written in the Lenard-Balescu 
form. We then extend this Lenard-Balescu part to < 3 so that its dynamically screened 
and Debye screened pieces may be treated separately. We thus have 

CT'''%Pa) = QT(Pa) + CT^'^iPa) • (B12) 

Here 

CZiPa) = Cil>'\p^) - A . j(-3)(pj (B13) 

corresponds to the result of our method of dimensional regularization: the first term 
^a6^^^(Pa) Boltzmann collision term for pure Coulomb scattering in u > 3 spatial di- 

mensions, with the scattering treated to all orders, while the second term involves ji5^^''(Pa), 
which is the number current of the Lenard-Balescu form (CI 5) with the dynamically screened, 
first Born approximation scattering amplitude given by Eq. (CI 7). The remainder reads 

^remain/ \ _ ^(;^>3)(1)/ \ , 9 t{'^<3)(1)/ \ i /^(''=3)(1) /'„ ^ (Ti^ A\ 

'-'ab \Pa) - -<^ab KPa) + ^ " -iuab KPo.) + <^oab IP«J • K^^V 

Here C'afe'^^'"'^'*(Pa) is the Boltzmann collision term for the pure Coulomb scattering in u > 3 
spatial dimensions in the first Born approximation, J^afe^'"'^''(Pa) is the Lenard-Balescu num- 
ber current for a Debye screened Coulomb potential in v <3 dimensions, and C'"^~^^'^^\pa) 
is the Boltzmann collision term for a Debye screened Coulomb potential in first Born ap- 
proximation in three spatial dimensions. Our final job is to show that the remainder (7^^mam 
vanishes, as we now shall do. 

If the difference of the Boltzmann collision terms that appear in Eq. (B14) is written as 
a single integral, then its integrand vanishes rapidly at high momentum transfer. The v = 3 
Boltzmann collision term for the Debye screened potential in first Born approximation may 
be extended to z/ > 3. To the appropriate leading order, this difference can then be expressed 
in the Lenard-Balescu form, 

-cr'""(pj + cn"«(pj = -A . jrs'(p.) , (B15) 

where 
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r(!^>3)(l) 
D—Cab 



Pa) 



k 7r(5 (va ■ k - Vfo ■ k) 



GaGb 


2 




1 


k?' + 




A;2 





k- 



k- 



d 



fa{Va)fb{Pb) ■ 



(B16) 



5P6 dPa. 

In the method of dimensional continuation, the pure Coulomb piece that appears here van- 
ishes. This is because the k integral carries the dimensions of L^~^ and there is no length L 
present to carry this dimension. The remaining term involving the Debye screened Coulomb 
potential defines, except for poles at integer dimensions, an analytic function for arbitrary 
dimensions that is identical in form with the current jj^aff'*'^^'' in the remainder (B14). Hence 



/^remain \ 



dPa 



r(i'<3)(l), 
^Dob ' 



Paj - 



dPa 



j(^>3)(l)/ N 



0. 



(B17) 



APPENDIX C: FOKKER-PLANCK AND LENARD-BALESCU 
LIMITS FROM BOLTZMANN EQUATION 

In general, if the squared scattering amplitude in the collision integral (Bl) decreases 
sufficiently rapidly at large momenta, or if the dimension v is sufficiently small such that the 
phase-space volume at high energies becomes small, then the collision integral (Bl) may be 
replaced by an equation of the Lenard-Balescu form. 

The Lenard-Balescu equation is a classical equation. The mechanical momentum transfer 
q and classical wave number k have the familiar relation 

q=/ik. (CI) 

It is the wave number k that is the significant variable in the scattering amplitude, while 
the momentum q appears in kinematical and phase space factors. The classical limit is the 
limit — > with k fixed. Thus we take the limit q — > in kinematical factors in which 
the momentum role is emphasized. This limit of small q in the kinematical factors in the 
Boltzmann equation produces the Fokker-Planck equation, to whose derivation we now turn. 
First we change variables for the 'spectator particle' h by writing 

1 1 
Pb = Pb + 2 *1 ' Pb = Pb - 2 *1 ■ (*^2) 

Removing the momentum-conserving delta function by the integration now presents the 
collision term (Bl) as 
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Cab{Pa) 



\Tf {2nh)6 



Pb ■ q 

rUb 



/fe Pb - fa{Pa + q) - /b P6 + /a(Pa 



(C3) 



We make an expansion in the momentum transfer q when it appears together with or in 
comparison with the momenta or p;,. As we shall soon see, the leading terms, the only 
terms that wc shall retain, are quadratic in q. 

Expanding to second order in q, and performing some rearrangement to simplify the 
result, produces 



1 



/b Pb - fa{Pa + q) - /b Pb + 7:q fa{Pa) 



1 



2 dpa 



d 



q- 



q- 



d 



fb{Pb)fa{Pa) ■ 



(C4) 



dpa dpb. 

Since this factor in the collision integral (C3) starts out linearly in q and since we shall work 
only to second order in q, the remaining terms need only be expanded to first order in q. 
Thus we may write 



Pa-q ^ q 



2m„ 



Pb • q 

nib 



1 1 d 
2 dpa 



Pa • q Pb • q 



m„ 



nib 



(C5) 



For a Galilean invariant theory, the scattering amplitude depends only upon the squared 
momentum transfer q^ and the (relative) energy in the center-of-mass W. By virtue of the 
conservation of energy [which is enforced by the delta function that remains in Eq. (C3)] 
this energy may be expressed in terms of either initial or final state variables. 



W = -mab (v„ - Vbf = - mab {y'a - ^'bf , (C6) 

where niab is the reduced mass and the velocities have the usual, generic, form v = p/m. 
Since p^ is the integration variable, we must write 

2 



W 



1 



rriab Va - Vb - 



2nib 



q 



(C7) 



Again, we need only consider the corresponding correction to linear order in q. According 
to the chain rule, this q-dcpcndence of W entails the derivative of |Tp with respect to W 
times 



q 



_d_ 

dVr 



W 



rriab 



(q-v„-q-V6) + C»(q)' 



(C8) 



2nib dva '^rub 
To leading order, this term involving the derivative of |Tp does not, in fact, contribute since 
the terms above multiply (5 (q ■ — q ■ Vb) and thus give a null result. Thus we may consider 
the scattering amplitude to simply be a function of 
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W^^ruab K-Vft)' . (C9) 

As we have remarked above, the squared wave number is the relevant variable for the 
case that we are now considering rather than the momentum transfer — h'^'k^. The usual 
form of the Boltzmann equation involves a squared scattering amplitude that is Galilean 
invariant and thus only a function of W and k^. However, we need to generalize this a 
little to take into account the plasma screening corrections that are included in the Lenard- 
Balescu limit of the Boltzmann equation. The background plasma specifies a rest-frame 
coordinate system, and so non-Galilean invariant variables may now also appear in the 
scattering amplitude. Since the system remains rotationally invariant, the only remaining 
variables involve the kinetic energies of the reacting particles. As we shall shortly see, the 
only relevant combination is the energy difference of the initial and final states of a particle 
as measured in the plasma rest frame. Since the total energy is conserved in the collision 
integral (C3), there is only one energy difference 

AE=lma {v'l -vl)=l m, [^l - v'l) . (CIO) 

Since 

A£^ = ^K + v„)-q, (Cll) 

this energy difference naturally defines a classical frequency, a frequency determined only by 
classical quantities, 

a; = A£;/n=^K + v„).k, (C12) 

and, to the order that concerns us, the squared scattering amplitude may be expressed*^^ as 

|T(iy,k2,K + q/2m„)-k)| ~ 

Since the expansion terms that are linear in q involve a complete integrand that is odd, they 
do not contribute, and the collision term (C3) now reduces to 

CabiPa) = • JabiPa) , (C14) 

in which 



1 1 d 
2 dpa 



|T(iy,k2,v„-k 



(C13) 



■^^As in the discussion about Eq. (C8), the derivative acting upon W gives no contribution to our 
order. 
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Jab(Pa) 



d 



k \hT{W, k\ Wa ■ k)r n6 (v„ • k - Vfe ■ k) 



k- 



dpb 



k- 



d 

dPa 



faiPa)MPb) ■ 



Here we have removed the overhne from the spectator momentum variable, Pb 
correspondingly written W in place of W. 

A trivial algebraic rearrangement of Eq's. (C14) and (C15) expresses the result as 



(C15) 
Pb, and 



Cab{Pa) = 



d 

dpi 



B[^{Pa) 



d 



KbiPa) 



f{Pa) , 



(C16) 



which is the Fokker-Planck form. 

The Lenard-Balescu equation entails only leading-order scattering, but fully dynamically 
screened. In this approximation. 



\hT\ 



GaGb 



Fe(k2,v„-k) 



(C17) 



and the reduced form (C15) becomes identical with the previous Lenard-Balescu coUision 
integral, Eq. (7.1). Note that the frequency dependence that appears here in the dielectric 
function, a; = Va • k is just the AE/h discussed above or, equivalently, the AE defined in 
Eq. (B6). Thus the scattering amplitude (C17) is just the first part of Eq. (BIO). 

As we have often remarked, the Lenard-Balescu equation is valid only for dimensions 
less than three, ly < 3. This reduced dimensionality is necessary for the convergence of 
the wave number integration at large k. In three dimensions, the dynamically screened 
Coulomb scattering amplitude does not vanish sufficiently rapidly at large q = /ik so as to 
permit the h expansion that we have made. If, however, the squared Debye Born amplitude 
is subtracted as in Eq. (BIO), then the Lenard-Balescu reduction may be made directly in 
three dimensions. Indeed, as we have repeatedly emphasized, this reduction must be made to 
consistently compute only the leading terms. Since the only kinematical variables available 
are momentum variables, the first correction, which is of order fi^, must appear in the form 
of a squared length = /i^/p^. A dimensionless ratio can only be obtained by multiplication 
with the Debye wave number. Hence the first correction to the leading classical limit that 
we have just derived is of order k^. This is formally of order of the plasma density 
n relative to the leading order result, a correction that is beyond the order to which we 
compute. Having said all this, one can now continue the Lenard-Balescu like equation in 



'Up to an omnipresent, omnivorous, logarithm. 
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three dimensions for the subtracted squared amphtude (BIO) to < 3. The two parts of the 
equation can then be separately treated as done in the previous Appendix B. 

Lifshitz and Pitaevskii [26] purport to derive the Fokker-Planck equation from the Boltz- 
mann equation for a dilute system of very heavy particles moving in a gas of light particles. 
Since our work might be confused as having some relationship to theirs, we briefly review 
it here. For the sake of completeness, we shall show that their work requires further ap- 
proximations to become internally consistent. This we do by providing an explicit example. 
But before starting out to do this, we should again emphasize that our reduction of the 
Boltzmann equation starts from the assumption that the scattering is restricted to small 
momentum transfers because of the dynamics of the scattering cross section. Lifshitz and 
Pitaevskii, on the other hand, assume that the momentum transfer is small because of the 
kinematics of a heavy particle moving in a light gas. 

First we transcribe the description of Lifshitz and Pitaevskii into our notation. They 
write the collision integral in the form 

CabiPa) = J d^qiWabiPa + Q, Q) fa{Pa + Q) " WabiPa, Q) /a(Pa)} • (C18) 

The "projectile" particle a is assumed to be heavy, and the "gas" particle b is assumed to 
be light. Comparing this structure with our standard form (Bl) and performing one of the 
momentum integrals trivially using the momentum-conserving delta function gives 

^abiPa) = J J^fM \T{<b,^')\" (27r;i)5 (^^q^ -q- v„,) . (C19) 

Here, as before, Vab = — is the relative velocity with = Po/"^a, '^b — P&/"^6, and 
niab is the reduced mass. Lifshitz and Pitaevskii expand Eq. (CIS) in powers of q when this 
momentum transfer appears added to the heavy particle momentum Pa and retain terms up 
to second order. This formal expansion gives the approximate collision term 

CabiPa) = ^ S^A^Pa) faiPa) + ^ ^^^(Pa) /a(Pa)} , (C20) 



where 



AUiPa) = J d'-q. WabiPa, q) , (C21) 



and 



^ir(Pa) -\ Jd'cL q^ WabiPa, q) ■ (C22) 
Lifshitz and Pitaevskii go on to set 
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AL{Va) = AIM + ^^^p ■ (C23) 

It is easy to see that the approximate coUision term (C20) vanishes when the heavy particle 
species a is in thermal equilibrium with the light gas particles b if 

AUPa)^PBlt{pa)vT. (C24) 

As we shall see, the satisfaction of this constraint requires further approximation. In this 
sense, the work of Lifshitz and Pitaevskii is misleading. 

The easiest and clearest way to demonstrate this is to consider an explicit example. 
We examine the case in which the scattering amplitude T is a constant and the spatial 
dimensionality is taken to be two, u — 2. These restrictions lead to trivial integrals and 
it is very easy to explicitly evaluate the expressions (C21) and (C22) using Eq. (C18) with 
the distribution of the light gas particles fb{Pb) of Maxwell-Boltzmann form at temperature 
T= The results are 

AUPa)^vlnb (Cmlb) , (C25) 

and 

BiTiPa) ={^TS^- + lm^b (vl V- + \ vl 5'-) } Ub {C m^,) , (C26) 

where Ub is the particle number density of the light gas particles, and C is a constant involving 
|Tp and geometrical factors such as 27r. These explicit results obviously violate the thermal 
equilibrium constraint (C24). A consistent result requires the further approximation that 
the "projectile" mass rria is much greater than the "gas" mass so that one can replace 
^ah/f^b 1- With this requirement obeyed, and assuming that the "projectile" speed is 
not too great in the sense that the kinetic energy mav'^/2 is not much greater than the 
temperature T, the second set of terms in Eq. (C26) may be neglected, and we see that the 
thermal equilibrium constraint (C24) is now obeyed. 

Be all this as it may, we must emphasize again that the Lifshitz-Pitaevskii treatment 
is based on a kinematical restriction which is very different than the dynamical condition 
that we impose in order to reduce the Boltzmann equation to a Fokker-Planck equation and 
ultimately obtain the Lenard-Balescu limit. 

APPENDIX D: THE CLASSICAL LIMIT 

Expression (10.13) for the cross section averaged momentum transfer is simply related to 
the classical limit. To see how this goes, we write Eq. (10.13) as 
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y"(icTg^ = 277^2 + 2 [1- COS 2 -(5;)] . (Dl) 

In the classical limit, large Z-values dominate. Thus, with the identification J = {l-\-l/2)h, the 
sum may be replaced by an integral, ti? 5^(^ + 1) — > / JdJ, and the phase shifts approximated 
by their WKB evaluation, 

= + dr yim{E -V)- Pjr''-'^ -pr^j , (D2) 

where E = /2m and is the turning point, the radial coordinate where the square root 
vanishes. The classical scattering angle 9 may be found in the usual way: The conservation 
of the angular momentum J = mr'^dO/dt is used to replace the time increment dt in the 
energy conservation equation by d,t = mr^d9/J. This gives the usual trajectory equation for 
dr/d9 which may be integrated to evaluate the classical scattering angle as 

/•oo T I 

^ = 7r-2 / dr— ^ (D3) 

Thus we find that, in the classical limit, 

e = 2h^ = 2 - 5i\ , (D4) 

so that in this limit 

jdaq^^27rj^ JdJ 2 {I - cosO) . (D5) 

The classical impact parameter h is defined such that J = p6, while = 2p^{l — cos9). 
Hence we have found that the classical limit of quantum mechanics yields 



oo 



daq^ = 27r / bdbq\ (D6) 
which is indeed the classical result. 
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